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Abstract 

We construct the effective Hamiltonian which governs the renormalization group flow of the 
gluon distribution with increasing energy and in the leading logarithmic approximation. This 
Hamiltonian defines a two-dimensional field theory which involves two types of Wilson lines: 
longitudinal Wilson lines which describe gluon recombination (or merging) and temporal Wilson 
lines which account for gluon bremsstrahlung (or splitting). The Hamiltonian is self-dual, i.e., 
it is invariant under the exchange of the two types of Wilson lines. In the high density regime 
where one can neglect gluon number fluctuations, the general Hamiltonian reduces to that for 
the JIMWLK evolution. In the dilute regime where gluon recombination becomes unimportant, 
it reduces to the dual partner of the JIMWLK Hamiltonian, which describes bremsstrahlung. 



Membre du Centre National de la Recherche Scientifique (CNRS), France. 
Preprint submitted to Elsevier Science 2 February 2008 



1 Introduction 



The construction of an effective field theory for scattering in QCD at high energy is a 
difficult, longstanding, problem, which has been at the heart of theoretical developments 
for more than twenty years. In the 'leading logarithmic approximation' (LLA) which is 
expected to control the high-energy limit, the effective action should include the quantum 
effects enhanced by the large energy logarithm r = In s/Q 2 (the 'rapidity'; as usual, s 
denotes the invariant total energy squared, and Q is the typical momentum transfer, or 
a particle mass), that is, it must resum radiative corrections of order (a s \ias/Q 2 ) n to all 
orders. If this is the only requirement, then the corresponding resummation is performed 
by the Balitsky-Fadin-Kuraev-Lipatov (BFKL) equation [I], which is however well known 
to violate unitarity in the high energy limit: the corresponding solution for the scattering 
amplitudes or for the gluon distribution grows like a power of the energy, in violation 
of the Froissart bound. According to the modern understanding, the BFKL equation 
should govern only the pre-asymptotic evolution, at intermediate values of the energy. 
The complete effective action valid in the high energy limit should unitarize the BFKL 
'pomeron' via the inclusion of multiple scattering. Alternatively, and equivalently, if the 
effective theory is written for the wavefunctions of the hadronic systems which participate 
in the collision (the 'target' and the 'projectile'), the effective action must include the 
non-linear effects, like gluon recombination, responsible for gluon saturation [2-5]. 

The first attempt towards constructing such an action has been given in Ref. [6], where the 
main arguments rely on kinematical and symmetry considerations. A lasting conclusion 
of this analysis is that the effective action should describe a two-dimensional field theory 
living in the plane transverse to the collision axis, and that the basic degrees of freedom 
should be Wilson line — path-ordered eikonal phases describing multiple scattering at 
high energy. The effective theory proposed in Ref. [6] involves two types of such Wilson 
lines, each of them describing the scattering of one of the participants in the collision off 
the color field created by the other participant. Subsequently, Lipatov and collaborators 
[7, 8] proposed an effective theory formulated in terms of reggeized gluons which contains 
the results in Ref. [6] as a special limit. More recently, Balitsky [9] has relied on a new 
factorization scheme (in rapidity) to construct an effective action in terms of Wilson lines, 
which extends the original construction in Ref. [6]. Some of his results will also emerge 
from our subsequent calculations. The direct comparison between the effective theory by 
Lipatov et al [7, 8] and that by Balitsky (or, more generally, any effective theory built in 
terms of Wilson lines, like the one to be developed below in this paper) is hindered by 
the fact that the general relation between the reggeized gluons and the usual gluon fields 
(and hence the Wilson lines) is not known. 

However, all the approaches mentioned above have merely considered the effective action 
corresponding to a given layer in rapidity : The action is obtained by integrating out gluon 
fluctuations with rapidity r within some intermediate range To < r < Ti, to leading order 
in ql s {ti — r ) and under the assumption that a s (ji — r ) <C I, in order for perturbation 
theory to apply. (The rapidity of a parton fluctuation is defined, as usual, as r = In l/x, 



2 



where x is the longitudinal momentum fraction carried by that parton.) From the point 
of view of the leading logarithmic approximation, this represents a single evolution step 
(in the presence of multiple scattering, though), but in order to study the high energy 
behavior one really needs to iterate this step. To that aim, one should be able to promote 
the effective action into a (field theoretical) renormalization group Hamiltonian which 
should describe the evolution of the effective theory with increasing r. In turn, this requires 
the proper identification of the relevant degrees of freedom and of the way that these are 
changed by the evolution in each rapidity step. 

Alternatively, the renormalization group (RG) flow can be formulated as an hierarchy 
of evolution equations for the relevant gluon correlation functions, which should reduce 
to the BFKL equation in the intermediate-energy regime where unitarity corrections are 
unimportant. Such an hierarchy has been derived for the first time by Balitsky [10]. This 
is an hierarchy for the evolution of Wilson-line operators which physically describe the 
scattering between a projectile made of elementary color charges (each of them identified 
by its own Wilson line) and the strong 'background' field of a target which is not evolving 
(like a large nucleus). The evolution consists in gluon splitting within the projectile wave- 
function, followed by the multiple scattering between the products of this splitting and 
the background field. With increasing r, the gluon density in the projectile wavefunction 
is rapidly increasing, yet this formalism does not allow for non-linear effects like gluon 
recombination, which on physical grounds are expected to tame this growth. 

The interpretation of Balitsky equations as an RG flow has been clarified by Weigert [1 1] , 
who noticed that the hierarchy in Ref. [10] can be compactly reformulated as a single, 
functional, evolution equation for a classical probability distribution for the Wilson lines. 
The functional differential operator in this equation plays the role of a Hamiltonian; it 
involves only one type of Wilson line (since the scattering is asymmetric: the target is 
dense and produces a strong color field, while the projectile is dilute), together with func- 
tional (Lie) derivatives with respect to these Wilson lines. A posteriori, one can recognize 
Weigert's Hamiltonian as a 'quantized' version of the effective action proposed by Balitsky 
for asymmetric scattering [10], where by 'quantization' we mean that the weak fields asso- 
ciated in the effective action with the projectile are replaced by functional Lie derivatives 
with respect to the Wilson lines built with the field of the target. 

Independently, a different RG approach aiming at the description of the gluon distribu- 
tion of a highly evolved hadronic system has been developed in Refs. [5, 12-15] and has 
eventually crystalized into an effective theory for gluon correlations at small x known as 
the color glass condensate (CGC) [15] (see also the review papers [16]). Following the 
general idea of the separation of scales in rapidity, this approach focuses on the evolution 
of a single hadron wavefunction (say, the target), with the purpose of integrating out 
the 'fast' quantum fluctuations — those with large longitudinal momenta — in layers of 
rapidity and thus obtain an effective theory for the correlations of the 'soft' (i.e., small-a;) 
gluons, as probed in a high energy collision. The Wilson lines appear naturally also in 
this approach: They describe the multiple scattering between the quantum gluons to be 
integrated out in the next evolution step and the strong color fields generated by 'color 
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sources' representing the fast gluons which have been integrated out in the previous steps. 

The central result in the CGC approach is the Jalilian-Marian-Iancu-McLerran-Weigert- 
Leonidov-Kovner (JIMWLK) equation [14, 15], a functional RG equation which turns out 
to be equivalent [15, 17] with the equation deduced by Weigert [11] from the Balitsky hi- 
erarchy [10]. This equivalence reflects the fact that both approaches include the same 
physical non-linear phenomena, which are treated either as recombination effects in the 
target wavefunction (in the CGC approach) or as gluon splitting in the projectile wave- 
function followed by multiple scattering off the target (in the approach by Balitsky). 
However, both approaches miss the other side of the problem, namely the saturation ef- 
fects in the evolution of the projectile (for the Balitsky equations) and, respectively, the 
correlations associated with gluon splitting in the target wavefunction (in the case of the 
JIMWLK equation). Because of that, both approaches miss the 'pomeron loops', as first 
recognized in Ref. [18]. These are the loops which open up with vertices for gluon split- 
ting (so like 2 — > n with n > 3) in the dilute regime and close through vertices for gluon 
merging (e.g., n — > 2) in the high density regime. (The denomination 'pomeron loops' is 
strictly correct only in the limit where the number of colors N c is large and the t-channel 
gluons can be pairwise combined into BFKL pomerons, but here we use this expression 
as a convenient short-cut also for the general case at arbitrary N c .) 

The gluon splitting and merging vertices in the high-energy evolution are generalizations 
of the 2-^2 BFKL vertex in the Lipatov Hamiltonian [19]. Such vertices have been explic- 
itly computed in perturbative QCD [20-24] (see [25] for a review) and appear implicitly 
as building blocks in the Balitsky and JIMWLK equations, from which they can be in 
principle extracted by expanding the Wilson lines in powers of the gauge fields. However, 
until recently the only formalism allowing for the simultaneous inclusion of both splittings 
and mergings in a high-energy scattering was Mueller's center-of-mass factorization of 
onium-onium scattering [26-28]. In this approach, the gluon splittings are included in 
the evolution of the wavefunctions of the two onia (described as evolving distributions of 
quark-antiquark 'color dipoles', as appropriate at large N c ), while the mergings are repre- 
sented by multiple scattering (the simultaneous scattering of several pairs of dipoles from 
the two onia). The numerical simulations of the onium-onium scattering by Salam [29] 
have demonstrated the importance of the correlations induced by gluon splitting, whose 
physical role has been fully understood only recently [30-32]: the gluon number fluctua- 
tions act as a seed for higher-point correlations in the dilute regime, and thus strongly 
influence the evolution towards high density. 

The first evolution equations in QCD to take into account both mergings and splittings 
have been proposed in Ref. [18] and shortly after improved in Refs. [33,34]. These equa- 
tions generalize Balitsky equations at large iV c by including the effects of gluon number 
fluctuations within the framework of the color dipole picture. By construction, these equa- 
tions describe dipole evolution (including dipole number fluctuations) in the dilute regime 
and JIMWLK-like gluon recombination in the high-density regime, and thus contain all 
the necessary ingredients to provide a complete picture of the evolution at large N c and 
for sufficiently high energies. The RG formulation of these equations is also explicitly 
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known [33,35]: this involves an effective Hamiltonian which describes a two-dimensional 
field theory for the dynamics of 'Pomerons' [35]. In this context, the 'Pomerons' repre- 
sent color singlet exchanges in the scattering between a dipole and a generic target. With 
increasing r, the Pomerons undergo BFKL evolution, they can dissociate (one Pomeron 
splitting into two) or recombine with each other (two Pomerons merging into one). This 
effective theory has a remarkable symmetry which reflects boost invariance: It is invari- 
ant under the duality transformations which exchange color fields with the functional 
derivative with respect to these fields (or, more precisely, with respect to their sources). 
Under these transformations, the splitting and mergings terms in the Hamiltonian get 
interchanged with each other, while the BFKL piece is self-dual. 

Whereas for large N c , the high-energy evolution in QCD is by now well understood 
[18,33-37], its generalization to finite N c remains an outstanding open problem, which is 
actively under investigation [38-40]. Recently, Kovner and Lublinsky [38] have constructed 
evolution equations for gluon correlations in the dilute regime valid for arbitrary N c , and 
then noticed [39] that the effective Hamiltonian generating these equations — to which 
we shall refer as the bremsstrahlung (BREM) Hamiltonian in what follows — is in fact 
dual to the JIMWLK Hamiltonian [11, 15]. Following this idea, they made the interesting 
suggestion [39] that the general (and yet unknown) evolution Hamiltonian, which describes 
both mergings and splittings, should be self-dual as a consequence of boost invariance. But 
in the absence of an explicit derivation of this Hamiltonian from the QCD Feynman rules, 
the hypotheses in Ref. [39] are difficult to test. Besides, by itself, the duality argument 
does not allow one to construct the general Hamiltonian for the QCD evolution, but only 
to relate its limiting expressions in the dense and dilute regimes. 

In fact, related work by the same authors [38, 40] shows that the extension of the JIMWLK 
equation to the dilute regime is afflicted with difficulties associated, in particular, with 
ambiguities due to the time-ordering of the correlations induced by quantum evolution. 
At this point, one should recall that previous formalisms like the CGC approach and 
the Balitsky equations, and also the effective action approaches in Refs. [6-9], have been 
specially tailored to deal with the high density, or strong field, regime, as relevant in 
the vicinity of the unitarity limit, but they need not be well suited for the description 
of subtle correlations induced by fluctuations in the dilute regime. Thus, in view of the 
current understanding, the very existence of a general evolution Hamiltonian which en- 
compasses all the desired phenomena is far from being clear, and the actual expression of 
this Hamiltonian, if any, is completely unknown. 

It is our main objective in this paper to investigate under which conditions the renormal- 
ization group approach developed in Refs. [14, 15] can be extended to account for gluon 
splitting in addition to gluon recombination, and thus promoted into a general theory for 
QCD evolution in the LLA. As we shall see, such an extension can be given indeed, but 
only within the limitations inherent to an RG analysis, which implies a coarse-graining as 
usual. Here, the 'coarse-graining' refers, of course, to rapidity, and thus to the longitudinal 
and temporal scales simultaneously. The effective theory can capture only those correla- 
tions that can be 'seen' by the quantum fluctuations which are integrated out in one step 
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of the evolution. As we shall briefly explain here (and then demonstrate at length in the 
main body of the paper), such correlations are fully encoded in two types of Wilson lines 

- path-ordered exponentials in the longitudinal and, respectively, temporal direction - 
which with increasing r evolve through gauge rotations at their end points. 

Indeed, one step of quantum evolution consists in integrating out 'semi-fast' gluon fluctu- 
ations in a rapidity layer r < r < T\ and in the presence of two types of background fields: 
the 'fast' fields created by the gluons with larger longitudinal momenta that have been 
integrated out in the previous steps, and the 'slow' fields radiated by the 'semi-fast' gluon 
itself. The multiple scattering off the 'fast' background fields accounts for gluon merging, 
whereas the simultaneous emission of several 'slow' fields implements gluon splitting (or 
bremsstrahlung) . The resolution scales of the semi-fast gluon are fixed by the light-cone 
kinematics together with the separation of scales: This gluon has a poor resolution both 
on the longitudinal scale characteristic for the fast fields and on the temporal scale char- 
acteristic for the slow ones. Thus, it can only measure an 'integrated' effect of these fields, 
and since this measure proceeds via multiple scattering, its result is encoded in Wilson 
lines. The supports of the longitudinal and temporal integrations implicit in the Wilson 
lines are evolving with r, but the background fields change only at the extremities of these 
supports, so the Wilson lines change only at their end points (where new correlations are 
induced by the quantum evolution). 

The effective Hamiltonian that we shall obtain is a functional of these two types of Wilson 
lines, and is moreover self-dual: that is, it is invariant under the exchange between lon- 
gitudinal and temporal Wilson lines. This invariance reflects a mirror symmetry between 
vertices for gluon splitting and, respectively, gluon merging, which is in turn related to 
the boost invariance of the evolution equations. In the limiting regimes where one of the 
background fields is weak (so that the corresponding Wilson lines can be expanded out 
in perturbation theory), the general Hamiltonian reduces to the expected limiting forms: 
(i) at high density, where the gluon number fluctuations become negligible, it reduces 
to the JIMWLK Hamiltonian [11, 15]; (ii) at low density, where gluon recombination 
can be ignored, it reduces to the bremsstrahlung Hamiltonian [39], which is dual to the 
JIMWLK Hamiltonian, in agreement with Ref. [39]. 

For each of the two limiting situations alluded to above (JIMWLK and BREM), the two- 
dimensional Hamiltonian structure is explicit and transparent: the canonical variables and 
the associated conjugate momenta are properly identified, and the corresponding Poisson 
brackets are explicitly known. Note that these Poisson brackets are non-trivial, in the 
sense that some of the variables — those associated with the weak fields and which play 
the role of Lie derivatives — are non-commutative: they do not commute with themselves, 
but rather obey the SU(iV c ) current algebra. The evolution equations are then obtained as 
the canonical equations of motion generated by the respective Hamiltonian via the Poisson 
brackets. For the JIMWLK Hamiltonian this procedure yields the Balitsky equations [10], 
as usual, whereas for the BREM Hamiltonian it leads to equations similar to those written 
down by Kovner and Lublinsky [38]. 
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But precisely because of the lack of commutativity alluded to above, we have not been 
able to extend the construction of the Poisson brackets to the general case, where both 
types of background fields are strong. The would-be straightforward construction, which 
consists in exponentiating the Lie derivatives, leads to non-local commutation relations, 
which reintroduce the longitudinal and temporal coordinates, and thus are unacceptable. 
The completion of the Hamiltonian structure in the general case is left for further studies. 

This paper will be organized as follows: In Sect. 2 we shall describe the scattering problem 
that we have in mind and its formulation in the CGC formalism. This will also give us 
the opportunity to introduce some notations and fix our conventions. In Sect. 3 we shall 
explain the general philosophy of the RG approach underlying the effective theory for the 
CGC. The discussion will be mostly qualitative, with the aim of anticipating and thus 
clarifying the physical interpretation of the subsequent, more technical, developments. 
The next two sections will provide the mathematical formulation of the RG: In Sect. 4 
we shall derive general formulae for the effective action encoding the correlations induced 
after one step in the quantum evolution. The action is manifestly gauge-invariant and is 
expressed in terms of the background field propagator for the semi-fast gluons, that we 
shall construct (within the approximations and for the physical regime of interest) in Sect. 
5. The self-duality will appear already in Sect. 4, as a symmetry relating two different, 
but equivalent, expressions for the effective action. The first application of the general 
formalism will be given in Sect. 6, where we shall provide a streamlined derivation of the 
JIMWLK Hamiltonian. This derivation is considerably shorter then the previous ones in 
the literature [14, 15], thus demonstrating the efficiency of the present formulation of the 
RG. Then, in Sect. 7 we turn to the general case and construct the effective action as a 
functional of the two types of Wilson lines. We shall derive four different, but equivalent, 
versions for this action. The self-duality is now manifest, as the invariance of the action 
under the exchange between longitudinal and temporal Wilson lines. We shall also check 
that the JIMWLK action is indeed obtained from the general action in the limit where 
fluctuations become negligible. Finally, in Sect. 8 we consider the limiting case of a dilute 
regime and thus deduce the Hamiltonian for bremsstrahlung. We construct the associated 
Poisson brackets and explain how to use them in order to derive evolution equations for 
observables like the scattering amplitudes. Sect. 9 contains our conclusions. 



2 The scattering problem in the CGC formalism 

Let us start with a description of the scattering problem that we have in mind. We shall 
consider a high energy scattering in which the target is a right-mover (it propagates in the 
positive z, or positive x + , direction), whereas the projectile is a left-mover (it propagates 
towards negative z, or positive x~). Note that, for the target, x + plays the role of time 
and x~ that of the longitudinal coordinate, while for the projectile the roles of x + and x~ 
are interchanged. Furthermore, we shall view the scattering in a Lorentz frame in which 
the target carries most of the total energy, so that its wavefunction is highly evolved - 
it contains many small-x gluons — and can be described as a color glass condensate (see 
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below). In that frame, the projectile is a relatively simple object, like a quark-antiquark 
pair in a colorless state (a 'color dipole'), or a collection of several such dipoles. 

Since the gluon density in the target is typically high, the projectile will undergo multiple 
scattering, which can be resummed in the eikonal approximation: A quark (or antiquark) 
which propagates through the target preserves a straightline trajectory but acquires a 
color precession; that is, its wavefunction gets multiplied by a Wilson line built with 
the projection of the color field of the target along the trajectory of the quark. For a 
left-moving quark with transverse coordinate x, this Wilson line reads : 

V\x) = P exp ji0 J dx~A+(x + ~ 0, x~, x) t a | , (2.1) 

where the t a 's are the generators of the SU(iV c ) algebra in the fundamental representation 
and the symbol P denotes the x~ ordering of the color matrices A+(x~)t a in the exponent, 
from right to left in increasing order of their x~ arguments. The integration runs formally 
over all the values of x~ , but in reality it is restricted to the longitudinal extent of the 
target, which is localized near x~ = because of Lorentz contraction. By the same 
argument, the projectile is localized near x + = 0, so it probes the target field A + at small 
x~^ , as also indicated on Eq. (2.1). 

In fact, the field A + is slowly varying in x + , since it is generated by fast moving color 
sources — the quarks and gluons inside the target — whose internal dynamics is slowed 
down by Lorentz time dilation. Thus, the target field can be treated as constant over the 
duration of the scattering. This motivates the color glass picture to be developed below, 
in which the scattering observables are first computed for a given configuration of the 
color fields in the target, and then averaged over the latter. 

Specifically, the S'-matrix for the scattering between a color dipole and the color field A + 
is obtained as the color average of the product of two Wilson lines: one for the quark and 
the other one for the antiquark (with transverse coordinates x and y, respectively) : 

S(x,y;A + ) = ±tr(ViV y ). (2.2) 

The physical S'-matrix for dipole-hadron scattering is finally obtained after averaging 
over all the configurations of the color fields A + in the target. Within the CGC effective 
theory, this average is computed as follows: 

(S{x,v)) T = ^-(tr(V;tvg) r ee J D[p\ W T [p]^-tr(V^V y ). (2.3) 

This formula makes it explicit that the field A + is created by 'color sources' (quantum 
fluctuations in the target wavefunction) with large longitudinal momenta k + , whose life- 
times are much larger than the collision time. In the computation of scattering amplitudes, 
these sources can be effectively replaced by a classical color current J£ = p a , where the 
charge density p a = p a (x~,x) is static, i.e., independent of x + , but random, as it corre- 
sponds to any of the possible configurations of the fast moving partons. The correlations 
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of p are encoded in the weight function W T [p] - a functional probability distribution 
which depends upon the target rapidity 2 r ~ Ins because of quantum evolution : 

The fast partons which represent the color sources have longitudinal momenta within the 
range P + ^> k + ^> A, where P + is the total momentum of the target and A is the scale at 
which the lifetime of a virtual excitation becomes of the order of the collision time. When 
increasing the energy s of the collision, P + increases like \/s, while A remains constant 
(since this is fixed by the properties of the projectile). Thus, the longitudinal phase-space 
P + k + ^> A available for quantum evolution is increased, leading to an enhanced gluon 
radiation, and thus to the generation of new color sources. The effects of this evolution 
will be discussed in detail and computed in the next sections. 

In particular, in the dilute regime at not so high energies, the target field is weak and 
Eq. (2.2) can be evaluated by expanding the exponentials in the Wilson lines. To lowest 
non-trivial order in this expansion one obtains the scattering amplitude in the two-gluon 
exchange (or single scattering) approximation: 

T(x, y) = 1 - S(x, y) ~ ^- (A+(x) - A+(y)) 2 , (2.4) 
where 

A+(x) = J dx~ A+(x + = 0, x~, x) (2.5) 

is the effective color field in the transverse plane, as obtained after integrating over the 
longitudinal profile of the target. Note that the ordering of the Wilson lines in x~ plays 
no role in this approximation, because of the symmetry of the color trace: tr(t a t b ) = 
l§ ab = tr (t h t a ). However, this ordering does matter in the strong field regime in which 
multiple scattering is important and the complete equation (2.2) must be used: successive 
collisions do not commute with each other (as they involve color matrices), so they resolve 
the longitudinal structure of the color field in the target. Thus in the construction of the 
effective theory one should also keep trace of the correlations of the color sources in x~ . 
But in the approximations of interest at high energy, these correlations turn out to be 
relatively simple and can always be embedded via Wilson lines. In particular, the effective 
theory and its evolution can be fully written in terms of Wilson lines [10, 11], in which 
case the longitudinal coordinate x~ is not explicit anymore. 



3 Quantum evolution: The general philosophy 

In this section we shall explain the general structure of the quantum evolution and de- 
scribe our strategy towards performing the calculation. This strategy relies on the renor- 

2 In the considered frame, r is essentially the same as the rapidity gap between the target and 
the projectile. 
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malization group (RG) in QCD at small-x, a technique which allows one to integrate out 
quantum fluctuations in the evolution with increasing energy by exploiting the separation 
of (longitudinal and temporal) scales specific to the high-energy problem. This technique 
generalizes the BFKL evolution [1] and has been originally developed in Refs. [14, 15] in 
relation with the JIMWLK equation. In what follows we shall further develop this tech- 
nique and apply it to the construction of more general evolution equations, from which the 
JIMWLK equation will emerge as a special limit. The discussion in this section is purely 
qualitative: we shall explain the structure of the quantum corrections and emphasize some 
general properties which will facilitate the physical interpretation of the calculations in 
the next sections. We shall thus anticipate from physical considerations some structural 
properties of the evolution equations that will be later confirmed by explicit calculations. 

As explained in the previous section, our objective is to construct an effective theory for 
gluon correlations in the target wavefunction at some soft longitudinal momentum scale 
A with A <C P + . For instance, to compute scattering amplitudes at the scale A, we need 3 
n-point functions like (A + (x)A + (y) ■ ■ -)\, where the subscript A means that the fields 
A + inside the brackets carry soft longitudinal momenta k + ~ A. At high energy, such 
soft fields are predominantly generated by relatively fast color sources, with longitudinal 
momenta p + such that P + 3> p + 3> A. Indeed, due to the infrared singularity of the 
bremsstrahlung spectrum oc dp + /p + , the contribution of these sources to the correlations 
at scale A will be enhanced by powers of the large energy logarithm ln(P + /A) ~ Ins, 
and thus will dominate over the genuine quantum fluctuations with k + ~ A. Therefore, 
in the leading logarithmic approximation (LLA) to which we shall restrict ourselves in 
what follows, the soft fields at scale A are entirely generated by fast color sources with 
p + ^> A and inherit the correlations of the latter. To the same approximation, the relevant 
'color sources' are gluons which are themselves radiated by other gluons with even larger 
longitudinal momenta. This separation of scales makes it natural to use an RG analysis 
in which the fast gluon fluctuations are integrated out in layers of p + and replaced by 
an effective color charge density p a (x~,x) which generates the same correlations as the 
quantum gluons at the soft scale A. 

Before we describe the quantum evolution, let us briefly review the solution to the classical 
equations of motion which rely the soft fields to p (this solution will be needed later). The 
relevant equations are the Yang-Mills equations with a source J£ = 5^ + p a , 



(D v F*">) a (x) = 5^ Pa (x), (3.1) 

where x = (x~, x) and D u = d v — igA a v T a with (T a ) bc = —if a bc- Because of the relatively 
simple structure of the current in the r.h.s., these equations can be solved explicitly (at 
least in specific gauges). Namely, it is consistent with Eq. (3.1) to search for a solution 
having the following properties [16] 



3 More precisely, we only need the correlations of the Wilson lines, like Eq. (2.3), but for the 
purpose of describing the quantum evolution it is more transparent to think in terms of color 
fields and their sources. 
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if = 0, A~ = 0, A+, A\ : independent of x + . (3.2) 

Then the only non-trivial field strength is F+\ After also imposing a gauge condition, 
the classical solution involves just one independent field degree of freedom. This becomes 
manifest in the Coulomb gauge V'i' = 0, where Eq. (3.1) reduces to 

-V 2 ± A+(x) = p a (x), (3.3) 

which is the light-cone analog of the Poisson equation. Note that we use a tilde to denote 
the classical source and field in the Coulomb gauge. In fact, the Coulomb gauge potential 
will appear so often in what follows that it becomes convenient to introduce a shorter 
notation for it: = a a . Eq. (3.3) is immediately solved as 



/cPy 1 f 
~T~ ln 7 ^r^Pa(x~,y) = / d 2 yA(x - y)p a (x~,y), (3.4) 
An (x — yyp A J 



where the infrared cutoff \i is necessary to invert the Laplacian operator in two dimensions, 
but it will eventually disappear from the physical results. The field strength in this gauge 
is obtained as F+ % = —d l a a . 

We shall also need later the classical solution in the light-cone (LC) gauge A^ = 0. This 
is constructed from the above solution in the Coulomb gauge via a gauge rotation 

A l {x-,x) = l -U{x-,x)d i U\x-,x), (3.5) 

with 

U ] {x~,x) = PexpS^igj X dz~ a a (z~ , a;)T a | . (3.6) 

where we use matrix notations appropriate for the adjoint representation (A 1 = A % a T a , 
etc). Note the boundary condition: A l a (x~,x) — > for x~ — > — oo. 

Eqs. (3.4), (3.5), and (3.6) together provide an explicit expression for the LC-gauge solu- 
tion A 1 in terms of the color source p in the Coulomb gauge. Note that the target average 
in equations like (2.3) can be computed by using the classical solution in any gauge, since 
the measure, the weight function, and the physical observables are all gauge invariant. 

We now return to the problem of the quantum evolution and start with the simplest 
case, that of the evolution of the 2-point function (2.4) in the dilute regime. In the LLA 
of interest, this evolution is described by the BFKL equation [1], and is schematically 
illustrated in Fig. I. 

Fig. l.a represents the 2-point function in the effective theory at scale A : the color 
fields A + are radiated from the classical charge density p and the average over p (with 
weight function Wa[p]) is represented by the upper blob. The one-step BFKL evolution 
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(a) (b) 

Fig. 1. The 2-point function in the effective theory (a) and its BFKL evolution in one step (b). 

of this 2-point function is represented in Fig. l.a: this rung diagram is a typical quantum 
correction which becomes important when we probe correlations at the lower scale 6A, 
with b <C 1. These correlations have two sources: the direct emission of a gluon with 
k + ~ bA from the classical source at scale A (cf. Fig. l.a) and the induced radiation 
from the semi-fast quantum gluons with intermediate momenta A ^> p + ^> bA (the rung 
in Fig. l.b). The quantum process gives a correction of order a s ln(l/6), and we shall 
assume that a s ln(l/6) <C 1 in order for perturbation theory to apply. The purpose of the 
renormalization group analysis is to replace the two contributions in Figs, l.a and b by a 
single classical contribution similar to Fig. l.a, but with a modified weight function W&aIv 9 ]- 
By computing the difference Wb\ — Wa ~ a s ln(l/6), one can establish the renormalization 
group equation (RGE) whose solution permits to iterate an arbitrary number of evolution 
steps. Introducing the rapidity r = ln(A/P + ), so that dr = ln(l/6) for one step of RG, it 
is clear that r plays the role of the evolution 'time'. 

It turns out that, quite generally, the RGE can be cast in Hamiltonian form, and can be 
applied either to the weight function 



^' iSp. 



Wr[p}, (3.7) 



or directly to the interesting observables, for example the dipole S'-matrix (2.3) 

8 



§p(S\p])r = -(H 



^' iSp 



Sip}) , (3.* 



where if is a Hermitian functional differential operator that we shall refer to as the 
'Hamiltonian'. By inspection of the diagrams in Fig. 1, it should be already clear what 
is the general structure of the BFKL Hamiltonian: When applied to a n-point function 
of the charge density p, this should annihilate two p's and replace them by two new ones 
which interact with each other with the exchange of a BFKL rung. This suggests that 

1 S S 

-^BFKL = - A^BFKL — PP~ , (3.9) 

2 op dp 

in schematic notations in which p should be understood as the two-dimensional color 
charge density, obtained as in Eq. (2.5) 
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p a (x) = dx p a (x ,X). 



(3.10) 



Indeed, there is no multiple scattering in the BFKL evolution, so the longitudinal structure 
can be trivially integrated out leaving a non-trivial dynamics in the transverse plane alone. 

Note two important features of Eq. (3.9) which are more general than the BFKL approxi- 
mation and will be recovered for the general equations below. First, the fact that the r.h.s. 
of Eq. (3.9) is a total derivative with respect to p is necessary to ensure probability con- 
servation: this guarantees that the normalization condition / D[p] W T [p] = 1 is preserved 
by the evolution. Second, note the specific ordering of the operators p and 5/5p (which do 
not commute with each other) : this is obtained after subtle cancellations between 'real' 
and 'virtual' corrections in the quantum calculation (see e.g. Refs. [11, 15]), but is in fact 
imposed by gauge symmetry [41]. Thus, in order to obtain the evolution Hamiltonian, it 
is in fact sufficient to compute the real correction (so like the real gluon emission in Fig. 
l.b) and then order the operators in such a way to ensure gauge symmetry. This is the 
strategy that we shall adopt in what follows. 

Let us turn now to the more interesting case where the target field is strong (this is 
the typical situation at high energy). This entails non-linear effects both in the classical 
equations of motion (3.1) and in the quantum evolution. Indeed, in this case, the semi- 
fast gluon in Fig. l.b propagates in the background of the strong classical field generated 
by the source p at scale A and undergoes multiple scattering off this field. Accordingly, 
the propagator of the semi-fast gluon must be computed to all orders in the background 
field. A typical diagram included in this resummation is shown in Fig. 2. a. Clearly, this 
describes gluon recombination : n gluons (with n > 2) are merging into two. Such merging 
processes are responsible for the saturation of the gluon distribution and the formation 
of a color glass condensate. 

In order to resum these processes to all orders, it is useful to notice that the kinematical 
conditions are satisfied for the use of the eikonal approximation . Indeed, as compared 
to the fast color sources responsible for the background field, the 'semi-fast gluons' are 
relatively slow. In the rest frame of the background field (which is the frame in which 
the eikonal approximation becomes most intuitive), the semi-fast gluons propagate in 
the negative z, or positive x~ , direction, so like the projectile in the scattering problem 
considered in the previous section. Therefore, the gluon propagator, and also the evolution 
Hamiltonian, will depend upon the background field via the same Wilson lines as in 
Eq. (2.1), except that these are now rewritten in the adjoint representation. 

We thus see that, in the presence of non-linear effects, the quantum evolution is able to 
discriminate the longitudinal structure of the color source, so like the scattering with an 
external projectile. This brings us to the issue of the x~ correlations generated by the 
quantum evolution. Because of the separation of scales in the problem, these correlations 
turn out to be quite trivial: they merely show that the color source p a {x~,x) extends 
towards larger values of x~ when increasing r [15]. Namely, by the uncertainty principle, 
the color source generated after integrating out the fast partons with momenta p + ^> A 
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(a) (b) (c) 

Fig. 2. One-step quantum evolution with gluon number changing vertices: ( a ) a n — > 2 merging 
process which contributes to the JIMWLK Hamiltonian; (b) a2^n splitting process as included 
in the BREM Hamiltonian; (c) a general m —* n process. 

is localized near x~ = within a distance Ax~ ~ 1/A. With the boundary conditions 
that we shall use in our calculation, and which are the same as in Ref. [15], the support 
of the source lies fully at positive values of x~ , namely at < x" < x~ with x~ ~ 1/A = 
(l/P + )e T . When a new layer of semi-fast gluons is integrated out, the additional color 
charge generated in this way is localized at 1/A ^ x~ ^ 1/&A, and thus has no overlap 
with the previously existing color source. Thus, the quantum evolution proceeds by adding 
new layers to the color source at larger values of x~, whereas the charge density in the 
inner layers — as generated in the previous steps — is not modified. This means that the 
correlations induced by the evolution are essentially local in x~ [17]. 

The previous discussion implies that, in a Wilson line like Eq. (2.1), the support of the 
integration is truly restricted to < x~ < x~ and, moreover, the net effect of increasing 
t by dr is an infinitesimal gauge rotation of the Wilson line at its large-x - end 

V}(x) = Pexp lw£ T dx~a(x-, a;) j — ► V} +dT (x) = J dxS* T (x- ^t^) ? 

(3.11) 

where ba a T has support in the outer layer at 1/A ^ x ^ 1/^A, but its detailed x 
dependence is irrelevant in the approximations of interest. Correspondingly, the functional 
derivatives in the RGE in the strong field regime are taken with respect to the color source 
(or field) in the outmost bin in x~ 

1 5 5 

-ffjiMWLK = - t — Xjimwlk[^, V*] - — . (3-12) 

2 oa T oa T 

This is the JIMWLK Hamiltonian whose first complete derivation (including the subtle 
issue of the r-dependence of the functional derivatives) has been given in Ref. [15]. Math- 
ematically, the functional derivatives which appear in Eq. (3.12) are really Lie derivatives 
with respect to the Wilson lines V and [17]. This is explicit in the independent deriva- 
tion of the JIMWLK Hamiltonian by Weigert [11] from the Balitsky equations [10]. 

Thus, in this strong field regime, the evolution Hamiltonian is naturally expressed in terms 
of Wilson lines and functional derivatives with respect to the latter, and not directly 
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in terms of the color source p a (x~,x). This is indeed the correct way to think about 
the high-energy evolution: as a renormalization group flow in a two-dimensional Hilbert 
space with a non-trivial geometry (namely, the group manifold spanned by the Wilson 
lines). On the other hand, given the separation of scales in the problem, one cannot 
address questions about the detailed longitudinal structure of the color source (like, e.g., 
computing correlations between the values of p a (x~,x) at different points x~) : only the 
correlations of the Wilson lines make sense. We shall recover this feature in more general 
situations below. 

Consider now the 2 — > n splitting diagram in Fig. 2.b, in which the n emerging gluons 
are all soft (i.e., they carry momenta k + ~ bA). Physically, this diagram describes gluon 
bremsstrahlung in the process of the BFKL evolution. Through such processes, n-point 
correlation functions with n > 2 get built from the 2-point function in the course of the 
quantum evolution. Gluon splitting has not been included in the analysis leading to the 
JIMWLK equation since formally this is a higher-order effect (as compared to the direct 
emission of n gluons from the classical source p) in the strong field regime. However, this 
argument implicitly assumes that the higher correlations for p have been already included 
in the weight function W r [p] . But when starting the evolution in the dilute regime (at low 
energy, or high transverse momenta), the higher-point correlations are originally absent, 
and they can only get built from the 2-point function. Thus, gluon splitting is in fact 
the dominant process for generating correlations in the dilute regime, and as such it has 
also a strong influence on the dynamics in the high-density regime, since the non-linear 
effects mix various correlations with each other, as obvious in Fig. 2. a. 

Focusing on the dilute regime for the time being, we would like to understand how to 
include bremsstrahlung in the BFKL evolution. To that aim, we notice that a diagram 
like Fig. 2.b can be interpreted as describing the propagation of the semi-fast gluon in the 
background of a soft color field (the radiated gluons). The background field is physically 
weak, as it corresponds to quantum radiation, but the associated effects will be resummed 
here to all orders (that is, we shall formally treat this field as being strong) since we want to 
generate all the 2 — > n processes for any n at once. Then the propagation of the semi-fast 
gluon can be again treated in the eikonal approximation. The longitudinal momentum of 
this gluon, comprised in the range A ^> p + ^> bA, is much larger than that of the radiated 
fields, so in the rest frame of the latter it will appear as a fast projectile moving in the 
positive x + direction. Accordingly, the semi-fast gluon couples to the minus component 
5A~ of the background field, via the temporal Wilson lines 



We have chosen to denote this background field as 5A~ to emphasize that it represents a 
quantum fluctuation at the soft scale bA, and not a classical field radiated by p (as previ- 
ously mentioned, the minus component vanishes for the classical solution; cf. Eq. (3.2)). 
Alternatively, in the scattering problem, 5A~ can be interpreted as the field radiated by 
the projectile, and which couples to the semi-fast gluon from the target. In the subsequent 
calculations, we shall generally use the first interpretation of 5A~ (as a quantum field), 




(3.13) 
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Fig. 3. The distribution of the various background Gelds in the x + — x~ plane. The overlap 
region around the tip of the light-cone is the interaction region relevant for quantum evolution 
(see Sect. 5). The contour surrounding the interaction region will be explained in Sect. 7. 



but both points of view will be useful in interpreting the final results. 

By definition, 8A~ has low k + and therefore large k~ (since this is a nearly on-shell gluon): 
k~ ~ Q 2 /2k + , with Q 2 a typical transverse momentum scale in the problem. Therefore, 
the field 5A~ is quasi-independent of x~ (more precisely, its respective variation can be 
neglected over the relatively small longitudinal extent Ax~ of the semi-fast gluon to 
which it couples) and is localized near x + = 0, within a small distance Ax + ~ l/k~. If 
5A~ is the field of the left-moving projectile, the same properties follow from the Lorentz 
contraction and the time dilation of the projectile. This explains why we have omitted 
the x~ variable in writing Eq. (3.13). 

The distribution of the various background fields in the longitudinal direction and in time 
for the effective theory at rapidity r is illustrated in Fig. 3. The color source p (and 
also the Coulomb field a or the field strength F +t ) has support at < x~ < x~ , with 
the upper limit x~ increasing with r. The soft field SA~ (and thus F~ l ) is localized at 
< x + < x+, where x+ decreases when increasing r. Indeed the product x+x~ ~ 1/Q 2 is 
constant under evolution. 

It should be clear by now that the result of resumming all splitting diagrams like Fig. 
2.b is quadratic in p and non-linear to all orders in SA~ , upon which it depends via the 
Wilson lines W and W^. It is furthermore intuitively clear, and will be more formally 
justified in the next section, that the effects of such diagrams on the evolution of the 
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n-point functions of the color charge p can be obtained with the replacement 



"•"^•"^W^r (3 ' 14) 

We thus anticipate that the evolution Hamiltonian describing bremsstrahlung in the dilute 
regime should have the following general structure: 

#BREM = ^PXBREM^W^p, (3.15) 

where p = p a (x) is integrated over x~ , so like in Eq. (3.10), since in this regime there is 
no multiple scattering in the x~ direction. What is, however, more subtle now is the x + - 
dependence of p, which was irrelevant in the BFKL and the JIMWLK approximations, but 
which becomes an issue in the presence of multiple scattering in the x + direction: With the 
prescription (3.14), the Wilson line (3.13) generates functional derivatives with respect 
to p a {x + , x) at all the points x + within the temporal extent Ax + of the fluctuating field 
5A~. Thus, the Hamiltonian (3.15) acts naturally on a Hilbert space in which the color 
charge density carries also x + dependence. This shows that, in order to include correlations 
associated with gluon splitting, one needs a generalization of the original CGC effective 
theory which allows for time dependence. We shall see later that such a generalization 
can be indeed given, but it cannot account for the detailed temporal correlations that are 
a priori encoded in a Hamiltonian like (3.15). This is a consequence of the separation of 
scales in the problem: the semi-fast gluon has small k~ relative to the background field 
SA~ , so it cannot resolve the structure of the latter in x + , but measures only its overall 
distribution, via the temporal Wilson line W. This is similar, but dual — in the sense 
of exchanging x~ by x + , A + by 5A~ , and thus V* by W — to the JIMWLK problem, 
in which the semi-fast gluon is unable to discriminate the detailed ^"-structure of the 
fast sources, but measures them only via the longitudinal Wilson line V'. Inspired by this 
analogy, one can already anticipate that, when interpreted in the sense of RG, the BREM 
Hamiltonian (3.15) describes the evolution of the correlation functions for W and W\ 
while the external factors of p act as Lie derivatives on these Wilson lines. This is what 
we shall find indeed in the subsequent analysis. 

The previous discussion suggests an interesting duality between the evolution in the dense 
and the dilute regimes, which in turn reflects a mirror symmetry between the n — > 2 and 
2 — > n vertices in Figs. 2. a and b, respectively. Specifically, the corresponding Hamilto- 
nians should correspond to each other via the following duality transformation : 

5 1 5 

x' < — > x + , a a (x-,x) < — > i — — — -, -— -< — > p a (x + ,x) . (3.16) 

dpa{x + ,x) % da a (x~, x) 

Anticipated by earlier works on a symmetric formalism of the high-energy scattering 
problem in QCD [6-8, 10], such a duality has been recently found to hold in the large-iV c 
approximation [35], where the only non-trivial vertices are the 1 — » 2 and 2 — * 1 triple- 
pomeron vertices, and the coordinates x + and x~ play no dynamical role. Furthermore, 
Kovner and Lublinsky have argued [39] that the duality between mergings and splittings 
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should hold exactly, i.e., beyond the large- N c approximation, as a consequence of a more 
general, self-duality, property of the general (unknown) evolution Hamiltonian. 

However, the authors of Ref. [39] have mistaken x + for x~ in the bremsstrahlung problem, 
and thus they failed to also recognize the x~ <-> x + duality. But after correcting for this 
misidentification, we shall see that the self-duality idea is essentially right. The general 
evolution Hamiltonian that we shall construct in the next sections turns out to be invariant 
under a generalized duality transformation which exchanges Wilson lines in x~ (so like 
V"t) with those in x + (like W). Moreover, the expressions for Hjimwlk and -^brem that we 
shall find are indeed dual each other in the sense of Eq. (3.16). Our explicit construction 
will further clarify the origin and the physical meaning of this duality. 

The general Hamiltonian alluded to above includes all the m — > n processes for arbitrary 
integers m > 2 and n > 2 (see Fig. 2.c). To resum such processes one needs the propagator 
of the semi-fast gluon in the presence of both types of background field: the field = a a 
generated by the fast partons {p + ^> A) and the field 5A~ associated with the soft quantum 
fluctuations (k + ~ 6A). This problem is complicated because the semi-fast gluons are at 
the same time slow relative to A\ and fast as compared to 5A~ , so in general the eikonal 
approximation cannot be applied. However, for the purpose of RG, the propagator is 
needed only for points x + and x~ near the origin of the light-cone, that is, in the region 
where the fields A+ and 5A~ overlap with each other (see Fig. 3). In this region, we shall 
be able to construct the propagator, and then the evolution Hamiltonian. The latter will 
be expressed in terms of Wilson lines (line-ordered phases in both x" and x + ) which run 
around the interaction region in Fig. 3 and represent the basic degrees of freedom. The 
general Hamiltonian turns out to be local in the transverse coordinates. The non-localities 
in the standard expressions for the JIMWLK [11, 15] and BREM [39] Hamiltonians (which 
here will be generated by expanding the general Hamiltonian in appropriate limits) are 
merely related to the non-locality of the relation (3.4) between the color field in the 
Coulomb gauge and its source. 



4 Quantum evolution: The effective action 



With this section, we start our explicit construction of the evolution Hamiltonian for 
QCD at high energy. This Hamiltonian must encompass processes like that in Fig. 2.c 
to all orders in the two types of background field : the ones above and those below the 
horizontal gluon line. The resummation of these diagrams can be formulated as a one-loop 
calculation in the quantum version of the effective theory for the CGC [14,15]. This is 
a theory in which the fast partons with longitudinal momenta p + ^> A are represented 
by a random color charge density p a {x~ ,x) with weight function Wa[p] while the slow 
gluons with momenta < A are still explicitly present as quantum fluctuations. Thus, this 
effective theory allows for the calculation of gluon correlations at any scale k + < A. 



18 



For instance, the 2-point correlation function is obtained as: 




D[p}W A [p} 



f K D5AA^{x)A v {y) exp(iS[A, p]) 
f A D5A exp{iS[A,p}) 



(4.1) 



where the symbol T denotes time-ordering and it is understood that the two fields inside 
the brackets 'live' both at the same soft scale, i.e., k± ~ k^ < A. Furthermore, the action 
S[A, p] describes the gluon dynamics in the presence of the 'external' (or 'background') 
color source p, and the total field A^ is the sum of the background field generated by p 
and the quantum fluctuations SA^ with momenta p + < A. As we shall see, the background 
field itself is influenced by the fluctuations. The action S[A, p] is such that it generates 
the classical field equations (3.1) in the limit where &4 M = 0, but its precise form will not 
be needed below. 

The one-step quantum evolution consists in integrating out the semi-fast fluctuations 
with momenta A ^> \p + \ ^> bA where b <C 1 but such that a s \n(l/b) <C 1. To that aim, 
we shall assume that the external fields in Eq. (4.1) carry momenta k + ~ bA and we shall 
compute the correlations induced at the soft scale bA after integrating out the semi-fast 
fields a M . That is, the total field within the action is now decomposed as 



where B^ is the background field generated by p, aF 1 represents the semi-fast fluctuations, 
and &4 M stands for the remaining fluctuations with k + < bA. Note that the total field at 
the scale bA (so like the external field in Eq. (4.1)) is simply N 1 = + 5A^. 

After separating the soft from the semi-fast modes in the functional integral, the quantum 
average in Eq. (4.1) becomes: 



where the last integration in Eq. (4.3) defines the effective action that we want to compute: 



This is an action for the soft fields &4 M which depends also upon p. From this action, the 
evolution Hamiltonian will be eventually obtained via the replacement (3.14) (or, rather, 
a generalization of it), that we shall shortly justify. 

So far, we have not specified the gauge in which the quantum theory is defined. With 
the separation of scales in k + being not gauge-invariant, this issue is indeed important. 
The partonic interpretation of gluons is meaningful in the light-cone quantization which 
is defined in the LC-gauge A + = 0, so this is the gauge that we shall adopt in what 
follows. Note that the separation of scales in k + is invariant under the residual gauge 
transformations in the LC-gauge, which are independent of x~ . In fact, the separation 
of scales can be also formulated in a gauge-invariant way, as we shall later explain, but 



A» = B» + a» + SA^, 



(4.2) 




(4.3) 




(4.4) 
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the use of the LC-gauge in the quantum calculation turns out to be really convenient. 
We should stress however that our calculation will be manifestly covariant with respect 
to the 'background' fields B^ and 5A^, so the final result for the effective action will be 
gauge-invariant . 

In the LC-gauge and for 5A^ = 0, the classical solution has only transverse components 
A 1 which are given by Eq. (3.5); these are time-independent and satisfy i™ = 0. But in 
order to compute diagrams like that in Fig. 2.c we need to preserve an explicit background 
field 5A~ (the component of the soft field SA^ which couples to the semi-fast gluon), and 
this will also modify the classical field problem, by introducing time-dependence. Namely, 
in the presence of SA~, the color source p will be subjected to color precession : 

p(x)=p a (x)T a — > J + (x+,x) = W x+ ,- 00 p(x)Wl + ^_ 00 , (4.5) 

meaning that the fast partons represented by p undergo eikonal scattering off the field 
5A~ . In the equation above, 

WV,-oo(a:) = P exp \ig dz + 5A-(z + , aj)| , (4.6) 

is the temporal Wilson line from z + — > — oo up to the actual time x + . Note that the 
current J M = 5^ + J + is covariantly conserved, as it should: 

D„J» = D~J + = (d~ -ig5A~)J + = 0. (4.7) 

We have anticipated here that B~ = 0, so that A" = 5 A" indeed. From Eq. (4.5), one 
sees that one can identify p with the asymptotic value of the current J + at x + — > — oo, 
and in what follows it will be useful to write: 

p-oo(x) = J + (x + — > -oo, a?) = p(x), Poo(x) = J + (x + — > oo,x). (4.8) 



The background field B^ 1 = 5^ l B l is defined as the solution to the Poisson equation with 
the current (4.5) : 



SS 



5A- 



= o =}► D V F V+ = J + = W 3& ,- 00 p- 00 (Z)Wl + _ oc , (4.9) 

B i +SA~ 



where, with a slight abuse of notations, the background field A^ = ^ % B % + 5^ 5A has 
been denoted as B l + 8A~ . 

A crucial observation for what follows is that 

F' + = -d + 5A~ ~ 0, (4.10) 

which follows from the fact that the field 5A~ contains only modes with very small longitu- 
dinal momenta and therefore it varies very slowly with x~ . More generally, the statement 
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that F + = represents the gauge-invariant expression of the separation of scales in the 
problem. With this property, Eq. (4.9) simplifies to 

D { F i+ = -D i d + B i = J+, (4.11) 

which will be shortly used to determine the field B l as a function of 5A~ and p. Notice 
the time dependence in the above equation : this is local in x + , so B % and J + have the 
same time dependence, which in turn comes from the color precession in Eq. (4.5). Since 
the soft gauge field 8A~ has support at small x + (cf. Fig. 3), the temporal Wilson line 
(4.6) looks like a smooth B-function in x + - - its variation with x + is concentrated in a 
small region above x + = 0. A similar x + -dependence then holds for both J + and B l . 

For the purpose of the one-loop calculation, the action S must be expanded to second 
order in the semi-fast fields a M . We obtain, with condensed notations, 



S = S[B* + 44", p] + H J+ l -a» G^[B l + 5A~,p] a" , (4.12) 

B l +&A- £ 



5A 

where 

5S 



5A l 
and 

5 2 S 



= D v F vi = DjF ji + D-F _i + D+F + \ (4.13) 



B l +SA- 



(4.14) 



As indicated in the equations above, only the SA component of the soft field 5A^ needs 
to be kept in the terms describing the coupling with a M . 

Using again the fact that d + 5A~ ~ 0, one can easily show that 

D_F- 1 ~ D + F + \ (4.15) 

a relation which lies at the origin of the duality property, as we shall later see. 

Consider now the Poisson equation (4.11) which determines B l . As compared to the 
corresponding equation for 5A~ = (the component p, — + of Eq. (3.1)), the present 
equation involves time dependence via J + . However, this dependence is local, so the 
solution to Eq. (4.11) is again of the form B l = (i/g)Ud % U^ (cf. Eq. (3.5)), except that 
now the Wilson lines depend also upon x + : 

U\x) = P exp ^9 J dz~a(x + , z~ , £c)| . (4.16) 

Here, a = B + is the background field in the Coulomb gauge, defined by the time- 
dependent generalization of Eq. (3.3) : 
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-V\a(x) = p(x) = U\x) J + (x)U(x), (4.17) 

which is local in both x~ and x + . The field B % constructed above is a two-dimensional 
pure gauge (F lj = 0), so Eq. (4.13) simplifies to 



SS 



5A* 



~ 2D + F +i ~ 2D_F~\ (4.18) 

Bi+SA- 



Performing the Gaussian integration over the semifast modes a in Eq. (4.4), we finally 
arrive at the following expression for the effective action: S e s — S + AS c s with 

AS«[6A-,p] = ~j (2D-F + XG i i(x,y) (2D~F +i ) b y , (4.19) 



xy 



or, equivalently (within the present approximations), 

AS cS [5A-,p] = -\\ (2D + F~r £ G^(x,y)(2D + F-%, (4.20) 



xy 



In these equations, iG l i are the transverse components of the time-ordered propagator of 
the semi-fast gluons: 

iG%(x,y) [B' + 5A-,p] = (T a£(x) a%(y)), (4.21) 

where the brackets denote the average over the quantum fields a M for fixed values of 
the background fields B l and 5A~ . This propagator can be obtained by inverting the 
differential operator in the r.h.s. of Eq. (4.14) in the LC gauge and with appropriate 
(Feynman) boundary conditions. 

Eqs. (4.19) and (4.20) are the fundamental formulae which will allow us to construct 
the evolution Hamiltonian in the next sections. As anticipated, these expressions are 
manifestly invariant under the gauge transformations of the background fields (with the 
propagator G u being covariant under such transformations). In what follows, we shall use 
this symmetry to replace the LC-gauge background fields = S^B 1 + <5 M ~&4~ by the 
Coulomb gauge ones A^ = 5 M+ a + 5^~5A , which are most directly related to the various 
Wilson lines. In the last equation, SA is the soft field in the Coulomb gauge: 



SA (x) = U\x) 



5A-(x + ,x) + -d- 
9 



U(x), (4.22) 



and depends on x (unlike its LC-gauge counterpart) via the corresponding dependence 
of the gauge rotation. That is, SA looks like a smooth 0-function in x~ . One can check 
on Eq. (4.22) that the condition F' + = (cf. Eq. (4.10)) remains satisfied, as it should. 

In fact, with the background fields in the Coulomb gauge, one has: 

D + = d + -iga, F +i = -&a, D~ = d~ - igSA~ , F" 1 = -d i bA . (4.23) 
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We shall later argue that, under the present approximations, the transverse propagator 

G&(x,y)[a,6A] = U\x) G^ b (x,y) [B\5A~] U(y) , (4.24) 
is invariant under the exchange of plus and minus components for coordinates and fields: 

x' < — > x + , a a (x) < — > SA a (x). (4.25) 

Therefore, under the transformations (4.25) the two expressions for the effective action, 
Eqs. (4.19) and (4.20), get interchanged. But since these expressions are equivalent with 
each other, we conclude that the effective action is invariant under the transformations 
(4.25). This symmetry is the most general form of the property that in the previous 
sections has been referred to as 'self-duality'. 

The physical meaning of this symmetry is easy to understand. The same one-loop calcu- 
lation resumming the diagrams illustrated in Fig. 2.c can be interpreted in two ways: (i) 
As the evolution of the color fields A + = a of the target with decreasing k + and in the 

presence of the fields SA of the projectile (this is the interpretation a priori privileged by 
our previous manipulations, and it corresponds to the expression (4.19) for the effective 
action), or (ii) as the evolution of the color fields 5A of the projectile with decreasing 
k~ and in the presence of the fields A + = a of the target (this interpretation corresponds 
to Eq. (4.20) and amounts to reading a diagram like Fig. 2.c upside down). By virtue of 
boost invariance, the two directions of evolution must be equivalent: the rapidity incre- 
ment dr can be given to either the target or the projectile, with identical results for the 
evolution of the scattering amplitudes, which are Lorentz-invariant. This equivalence is 
indeed manifest at the level of our results, as the invariance of the effective action under 
the change (4.25) of the direction of evolution. 

To conclude this section, let us justify the replacement (3.14) which will allow us to 
transform the effective action (a functional of p and SA~) into an evolution Hamiltonian 
involving p and 5/5p. This step is essential for promoting the above one-loop calculation 
into a systematic renormalization group (RG) procedure. The purpose of RG is to shift 
the newly induced correlations at the soft scale bA from the effective action AS c g to 
the weight function W T [p] for p. To that aim, it is preferable to work in the Coulomb 
background gauge, since scattering observables are directly related to A + = a, and not 
to B\ Consider then the evolution 

(Ta(x)a(y)) bA - (Ta(x)a(y)) A (4.26) 

of the 2-point function when decreasing k + from A to bA. We would like to associate 
this evolution with a change a — > a + 5a in the classical field ; through Eq. (4.17), 
this corresponds to a change p — > p + 5p in the charge density in the Coulomb gauge. 
Physically, Sp represents the color charge density of the semi-fast gluons. To reproduce 
the correlations encoded into AS e g, the induced source Sp must be a random variable (for 
a fixed distribution p of the color charge at scale A), with correlations: 
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(T5p(x 1 )5p(x 2 ) ■ ■ ■ Sp(x n ))' 



,-n+l 



^nyEli ■ ■ ■ i -En) — 



,-n+l rn 



M-(xi)M-(x 2 )...M-(a; n ) 



A-=0 

(4.27) 



where for more clarity we have temporarily denoted 5A simply as A . Indeed, the quan- 
tum calculation yields (for the example of the 2-point function, once again) 



(T5A + (x)5A + (y))~ = j^DSA 5A + (x)5A + (y) exp {i(S[5A} + AS eS [SA , p])}, (4.28) 

where the path-integral in the r.h.s. is now evaluated in the Coulomb gauge. To the 
accuracy of interest, this has to be computed in the saddle point approximation and to 
order a s ln(l/6) (recall that the induced action AS eS is itself of 0(a s ln(l/6))). This yields: 



(T5A + (x)5A + (y))~ ~ i J iG + -(x,u)iG + -( yi v) F 2 (u,v) , 



(4.29) 



which should be compared to the corresponding correlation induced, via the classical field 
equation (4.17), by a change Sp in the color source: 



(T5a(x)5a(y))' 



A(x - u)A(y - v) (T 8p{u) 5p(v))~. 



(4.30) 



In the static limit p =0, G + (x,u) = A(x — u), and the above correlations coincide 
with each other provided Eq. (4.27) is satisfied. 

Finally, it is clear that the Hamiltonian which, when acting on correlations of p (cf. 
Eq. (3.8)), generates the correlations in Eq. (4.27) is of the form (with dr = In 1/6) : 



-drH 



P, 1 



6_ 

Tp 



Y^—.\ r n (x 1 ,x 2 ,...,x n ) 

Ail 

= iAS Pff p, SA 



i n+1 5 r - 



Sp(x 1 )Sp(x 2 )...Sp(x n ) 



5 



5p 



(4.31) 



Thus, the evolution Hamiltonian is obtained from the effective action via the replacement 

5 



5A (x) 



6p(xY 



(4.32) 



which generalizes the previous prescription (3.14) by including the dependence upon x~ 
and by specifying that in the general, strong field, regime 4 , this replacement works at 
the level of Coulomb gauge quantities. 



4 Recall that Eq. (3.14) has been written for the dilute regime, in which the distinction between 
the LC-gauge and the Coulomb gauge for the background fields becomes irrelevant, since one 
can approximate U « « 1 in any gauge rotation like Eq. (4.22). 
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5 The background field propagator 



In this section we shall construct the background field propagator which is needed to 
compute the expressions (4.19) or (4.20) for the effective action. This is the propagator 
G l ^{x ) y) for a semi-fast, transverse, gluon in the quantum LC-gauge (a + = 0) and in 
the presence of a background field which can be taken either in the LC-gauge (A^ = 

8^ l B l + 5^~5A~), or in the Coulomb gauge (A^ = 5^ + a + 5^~5A ), and which involves 
two independent field degrees of freedom. The general expression of this propagator is 
complicated and we shall not attempt to derive it here. But for the present purposes, a 
crucial simplification comes from the fact that Eqs. (4.19)-(4.20) involve the propagator 
G l i(x, y) only at points x and y within the interaction region in Fig. 3. By the 'interaction 
region' we mean the diamond-shape area around the tip of the light-cone where the 
supports of p and 5A~ (or, equivalently, of a and 8A ) overlap with each other. The 
essential property of this region is that its extent both in longitudinal and in temporal 
directions is very small as compared to the respective resolution scales of the semi-fast 
gluon. Because of that, the propagator G^(x,y) is local in the transverse coordinates and 
has a relatively simple structure in the other coordinates, as we shall explain below. 

Let us first establish the support properties of the integrands in Eqs. (4.19) or (4.20). 
Consider Eq. (4.19), for definiteness, and look at the vertices there: the field F +l has the 
same support as the color source p at scale A, so it is localized in x~ within the region 
^ x~ < 1/A. Furthermore, its (covariant) derivative D~ F +l = (d~ — ig5A~)F +l involves 
two terms, which are both localized at small x + , within the support of 5A~. This is obvious 
for the second term, which is explicitly proportional to SA~ , but it is also true for the first 
term d~F +1 , since the time-dependence of F +l comes entirely from its interaction with 
5A~ , as discussed after Eq. (4.11). Recalling that SA~ has small k + ^ bA and thus large 
k~ ^ Q 2 /2bA (with Q 2 a typical transverse momentum scale), we deduce that D~ F +l has 
support at ^ x + ^ bA/Q 2 . To summarize: 

for 0<:r~<l/A and < x + < bA/Q 2 . (5.1) 

This is also the support of the vertex D + F~ l in Eq. (4.20). On the other hand, the 
semi-fast gluon has A ^> k + ^> bA, and thus typical resolution scales Ax~ ~ 1/bA and 
Ax + ~ A/Q 2 . Since, by assumption, b 1, we conclude that the semi-fast gluon cannot 
discriminate the internal structure of the interaction region in Fig. 3, as anticipated. 

This feature drastically simplifies the calculation of the semi-fast propagator. To better 
appreciate that, let us analyze first the corresponding free propagator 5 : 

CTS(p) = ^Go(p), G* -(p) = ^—G (p), G-«(p) = _£!_G (p), (5.2) 

where G (p) = l/(2p + p~ — p 2 + ie) is the same as the propagator of a free scalar field. 
5 The component Gq~ is not shown since it is irrelevant for the present purposes. 
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Consider the scalar propagator, and perform the Fourier transform to coordinate space: 



/ffff^' —ip^ (x — y J r 
Ji_ e [Q(p+)Q(x + - y + ) - Q(-p+)Q(y+ - x + ) 

strip 

where the integration over p~ has been computed via contour techniques and the in- 
tegration over p + is restricted to the strip A ^> \p + \ ^> bA. The interesting region is 
where both x~ — y~ and x + — y + are small compared to the typical variation scales, l/p + 
and l/p~ = 2p + /p 2 , respectively. Then we can neglect the longitudinal and temporal 
coordinates within the exponentials and perform the remaining integrations, to obtain: 

iG l j (x-y)~ 5 lj -i- In \ 5 {2) (x-y), (5.4) 

4:71 b 

where the logarithm has been generated by the restricted integration over p + : 



rrfp+0(p+) = _ r dp+ Q(-p+) = J_ 1 

7 2vr 2p+ 7 2vr 2p+ 4tt 6' 1 ' J 

strip strip 



The final result, Eq. (5.4), is the same as the Fourier transform of the simplified momentum- 
space propagator in which p 2 is neglected compared to p + p~ : 

strip 

This is intuitive since small values for x + and x~ correspond in momentum space to large 
values for p~ and p + , respectively. Similar approximations can be performed also in the 
presence of the background fields. 

Consider first the case where there is only one such a field, namely 5A~. In the LC- 
gauge a + = 0, insertions of 5A~ cannot couple the transverse components to any 
other component like G l ~\ thus, the calculation of G y [SA~] is tantamount to that of the 
respective scalar propagator. The general expression for G t: '(x,y)[5A^] valid for arbitrary 
points x and y can be inferred from previous results in the literature (see, e.g., Ref. [15]). 
Here, however, we only need its expression valid within the interaction region, and this 
can be easily computed by solving the following differential equation 

2dt[d~ -tg5A-(x)]G(x,y) =6^(x-y), (5.7) 

which has been obtained from the general Dyson equation after neglecting the transverse 
Laplacian relative to d + D~, so like in Eq. (5.6). Working in the p + representation 
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(since SA is independent of x ) and using D x W x + y +(x) = for x + > y + , where 

W x+y+ {x) = Pexpjic/^ dz + 5A~(z + ,x)^ , (5.8) 

one immediately finds that the solution to (5.7) with Feynman boundary conditions reads 

tG^(x + ,x,y + ,y,p + )[5A-] = S ij ^-S^(x - y)[<?>(p + )Q(x + - y + )W x+y+ (x) 

zp 



(5.9) 



-Q(-p + )Q(y + -x + )Wl +x+ {x) 

After also performing the restricted integration over p + , one finally obtains : 
iG^(x + ,x,y + ,y)[5A-} = 5* In \ 5 {2 \x - y)\Q(x + - y + )W x+y+ (x) 

Quo 1 

+ Q { y+_ X + )W l +x+{x) ^ (5.10) 

Consider similarly the case where the only background field is the one created by p, which 
is independent of x + . If the field is taken in the Coulomb gauge, where it has just a plus 
component B + = a, then the corresponding propagator is obtained by simply replacing 
x + — > x~ and 5A~ — > a (and therefore W x + y +(x) — > U x - y - (x)) in Eq. (5.10): 

i&i(x-, x, y~, y)[a] = 6* i- In I ^(z - y) [e(x" - y")^-,- (*) 

+ e( J /--a;-)C/ v - a: -(x)]. (5.11) 

Note that, a priori, the insertions of a can couple the transverse semi-fast gluons a 1 
to the temporal ones a~ , so in general the propagator (^[ck] will involve a mixture of 
all the components of the free propagator shown in Eq. (5.2). Still, for small x~ and y~~ 
(within the interaction region) the coupling between transverse and temporal components 
is negligible since the off-diagonal propagators like G l ~ are suppressed when p + ^> p\ 

If one now rotates the propagator (5.11) to the background LC-gauge, — > G^[B l \ 
with G^(x,y) = U(x)G i: '(x, y)W(y), cf. Eq. (4.24), then one discovers that within the 
present approximations G l ^[B l ] is the same as the free propagator (5.4) ! This can be 
easily understood: Unlike the Coulomb gauge field a which is singular at x~ = (on 
the resolution scale of the semi-fast gluons), the LC-gauge field B % is only discontinuous 
there, so the corresponding propagator G l ^[B l \ must be continuous at both x~ = and 
y~ = 0. Thus, for small separations x~ — y~ ~ 0, the propagator will be independent of 
x~ and y~, and therefore also independent of the field B l . 

Armed with this experience from the simpler cases, we now turn to the general situation 
where both types of background field are simultaneously present. We shall first take these 
fields in the (background) LC-gauge, where A^ = 5^B l + 5^~5A~. Then the propagator 
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G l ^[B l , 5A~] will be continuous at small x~ (and thus independent of B l ), but discontinu- 
ous at small x + (and thus dependent upon 5A~ via the corresponding Wilson line (5.8)). 
We conclude that the general propagator that we need for the present purposes is the 
same as the propagator (5.10) in the presence of SA~ alone. This is the propagator that 
we shall use for our most general calculations below. 

But it is still interesting (in particular, for the self-duality argument in Sect. 4) to also 
consider the general propagator in the background Coulomb gauge, where A 11 = S^a + 
5^~5A . This is obtained from Eq. (5.10) via the gauge-rotation (4.24), and it is found to 
involve products of temporal and longitudinal Wilson lines, of the type: 

W x+y+ (x-,x)Ul- y 4y + ,x) = Ul- y -(x + ,x)W x+y+ (y-,x), (5.12) 

where W x + y + is the temporal Wilson line in the Coulomb gauge: 

W x+y +(x~,x) = U jf (x~,x + ,x) W x+y +(x)U(x~,y + ,x) 

= Pexp\ig [ X dz + 5A~ (z + , x~ , x)\ , (5.13) 



with SA the corresponding field, cf. Eq. (4.22). The equality in Eq. (5.12) follows because, 
in the LC-gauge, W x + y + is independent of x~ , by Eq. (4.10). Note that Eq. (5.12) involves 
two possible paths joining the points (x~ , x + ) and (y~ ,y + ); both these paths give the same 
contribution because of the condition F H = 0, which reflects the separation of scales in 
the problem. But by the same condition, any other path joining these points will give an 
identical contribution, so the expressions in Eq. (5.12) can be equivalently replaced by 

Pexpj^y (dz + 6A~ + dz~a)(z + ,z~,x)^ , (5.14) 

where the path 7 going from (y~,y + ) to (x~,x + ) must lie inside the interaction region, 
but otherwise is arbitrary. With the Wilson line (5.14), the Coulomb gauge propagator 

G lj [a, SA ] is clearly invariant under the transformations (4.25), as anticipated in Sect. 4. 



6 The high density regime: JIMWLK Hamiltonian 



Before we address the general case in the next section, it is instructive to see how the 
formalism works in a case for which the final result is already known: the strong field/no 
gluon number fluctuation regime, where the RG evolution is governed by the JIMWLK 
Hamiltonian [11, 14, 15]. For more variety in the presentation, the path that we shall follow 
in this section to recover the JIMWLK Hamiltonian will be different from that to be used 
in the next section in relation with the general case. 

As it should be clear from a comparison between the diagrams in Figs. 2. a and c, and also 
from the general discussion in Sect. 3, the effective action corresponding to this regime 
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is obtained by retaining terms to quadratic order in 5A~ and to all orders in B l (or a) 
in the general expressions (4.19) and (4.20). Since the vertices D~F +l or D + F~ % in these 
expressions start already at linear order in 5A~ , it will be enough for the present purposes 
to consider the expression of the transverse propagator for 5A~ = 0. It turns out that the 
JIMWLK Hamiltonian is most efficiently obtained starting with Eq. (4.20) and working 
in the background Coulomb gauge. That is, we shall start with 

AScS = ~\l ( 2 D + F-%&(x,y)(2D + F-%, (6.1) 

xy 

where D + = d + — iga, F~ l = —d l 5A , and the propagator G tj is given by Eq. (5.11). 

The first step is to integrate by parts over the longitudinal coordinates x~ and y~. From 
Eq. (5.11), one can easily check that D+G tj (x,y) = 0. By using this, together with the 
identity (A and B are two generic color matrices) 

&*(AB) = (D»A)B + A(D^B). (6.2) 

one can immediately deduce that the integrand in Eq. (6.1) is a total derivative w.r.t. 
both x~ and y~ , so the result of the integration comes from the 'surface' terms at large 
(positive or negative) longitudinal coordinates. Specifically, 

AS cS = In 1 J [FJ(oo, x)F£(oo, x) + F^(-oo, x)F y+ \-oo, x) 

x+,y+,x 

- F- + l (oo, x)V\x)F; + %-oo, x) - F-^-oo, x)V(x)F; + %oo, x)] , (6.3) 
where (cf. Eq. (3.11)) : 

V\x) = Pexp ji0 J dx-a a (x-,x)T a } . (6.4) 

Note that the time arguments of the fields in Eq. (6.3) are shown as lower subscripts, a 
convention that we shall often use in what follows. As explained in Sect. 4, the difference 
between the fields F~ l (x~ = oo) and F~ l (x~ = — oo) comes from the interaction with 
the color source p, as encoded in the gauge rotations in Eq. (4.22). This interaction is 
localized at small x~ , so F~ % (x~) looks like a smooth 6-function in x~ . 

By using F~ % = —d l 5A and defining A~ (x) = J dx + 5A (x), we obtain 

AScff = lni / \d i A-{oo,x)d t A'{oo,x) + d i A^{-oo,x)d i A-(-oo,x) 
2n b J L 



X 



- &A-{po, x)V\x)d i A-{-oo, x) - &A-(-oo, x)V(x)d i A-(oo, x)] . (6.5) 

The evolution Hamiltonian can be now obtained via the replacement (4.32) (in four- 
dimensional coordinates !). By also using the Poisson equation (4.17) to reexpress S/Sp 
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in terms of 5 /5a, we finally arrive at 

-f^JIMWLK = ?2~T3 J K-xyz 
xyz 



+ 



5a a (oo, x) 5a a (oo, y) 5a a (—oo, x) 5a a (— oo, j/) 



<5a: a (oo, a?) a 5a b (—oo,y) 5a a (—oo,x) 5a b (oo,y) 
with the following transverse kernel: 



(6.6) 



IC(x ,y, z) = J — W( — w . (6.7) 

In writing Eq. (6.6) we have suppressed the x + dependence although, strictly speaking, 
the integrand does depend upon time, via the functional derivatives which should be 
interpreted as, e.g., 

- = f dx + - . (6.8) 

5a a (oo,x) J 5a a (x + ,oo,x) 

But this dependence is rather trivial since when the Hamiltonian acts on scattering ob- 
servables, like Eq. (2.1), the only contribution arises from the interaction time x + ~ 0. 

On the other hand, the x~ dependence of the Hamiltonian (6.6) is less trivial and deserves 
some comments: As mentioned before, 5A (x) looks like a smooth O-function in x~ , which 
varies within the range ^ x~ ^ x~ , with x~ ~ 1/A = (l/P + )e r . Therefore, in the 

identification (4.32), the asymptotic fields 5A (x~ — > ±oo) should be really understood 
as the values of the field at the end points Xq ~ and x~ of this finite interval. Thus, 
strictly speaking, 

5 5,5 5 

— — t r = — — -—- and — — t = — — (6.9) 

oa a {oo,x) oa^{x) oa a {—oo,x) oa^x) 

act as derivatives at the end points of the Wilson line in Eq. (3.11): 

By comparing these results and using the identity V ] ab T b = VT a V\ one can express the 
derivatives at Xq in terms of those at x~: 

V b [(x) = V ab {x) -J-r . (6.11) 



5a%{x) 5a b T {x) baK ' auy ' 5a b {x) 

Alternatively, and more generally, Eq. (6.11) can be obtained as the result of a gauge 
rotation, as we show now. Let us first introduce the gauge-covariant 'color current' : 

J-(x) = D v F v -{x) = DiF l ~ + D + F + - = D l F~\ (6.12) 
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where in writing the last equality we have used F + ~ = 0, cf. Eq. (4.10). (In a scattering 
problem, in which 5A~ would be the field created by a left-moving projectile, J~ would 
be the color current associated with the latter.) Under the gauge rotation (4.16) from the 
LC- to the Coulomb gauge, this current transforms like: 

J-{x) = U\x) J - (a;) U{x) , (6.13) 
a relation which for x~ = oo can be rewritten as: 

J-(x + ,oo,x) = U j ab (x + ,oo,x) J-(x + ,-oo,x) . (6.14) 

By also using F~ l = — d^A together with Eqs. (6.12), (4.32) and (4.17), one can succes- 
sively write 

dp a (x) da a (x) 

By inserting this representation for J~(x) in both the l.h.s. and the r.h.s. of Eq. (6.14), 
and recalling that W(x + , oo, x) ps V\x) in the present regime (small 8A~), we finally 
recover the relation (6.11), as anticipated. 

After using Eqs. (6.9) and Eq. (6.11), the JIMWLK Hamiltonian (6.6) takes its standard 
form: 



a JIMWLK 



(6.16) 



xyz 



5a b T (y 



as originally derived in Ref. [15]. Note that, although we have explicitly computed here 
only the 'real' quantum correction, i.e., the gluon exchange diagram in Fig. 2. a, the above 
Hamiltonian is complete as written, that is, it also contains the 'virtual' correction (asso- 
ciated with self-energy and vertex corrections), which is generated when commuting the 
functional derivatives in Eq. (6.16) through the Wilson lines in the kernel there [11, 15]. 
The proper ordering of the operators inside the Hamiltonian (which implicitly takes into 
account the virtual corrections) has naturally emerged from the previous manipulations 
because of the intimate connection between the ordering of the operators and the gauge 
symmetry : The operators in Eqs. (6.6) or (6.16) are line-ordered in x~ in such a way 
that the Hamiltonian be invariant under the gauge transformations which depend only 
upon x~ (the residual 6 gauge transformations in the Coulomb gauge). 

Let us conclude this section with some considerations about the Hamiltonian structure 
associated with the JIMWLK Hamiltonian (see also Refs. [11,15,17,41]). Such consider- 
ations may look unnecessarily formal at this level, but they will be useful for clarifying 
more general cases later on. So far, we have implicitly assumed that -£/jimwlk acts as a 



6 These are the transformations which preserve the structure of the classical field in the Coulomb 
gauge: = 5^ + a with a independent of x + ; see Ref. [41] for details. 
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(second-order) functional differential operator in the Hilbert space spanned by the ob- 
servables 0(V, V^) built with the Wilson lines. This is sufficient for the present purposes 
since, as we have seen in Sect. 2, all the observables pertinent to the scattering between 
the high-density right-moving target and a (relatively simple) left-moving projectile are 
built in terms of Wilson lines. Within this Hilbert space, the two functional derivatives in- 
troduced in Eq. (6.9) act, respectively, as left and right Lie derivatives 7 (i.e., as generators 
of the infinitesimal left and right gauge rotations). Specifically, if one writes 



i 8a^(x) ' 



= V ab (x) J b L {x) , 



(6.17) 



then one has, for instance (with the Wilson lines taken in some arbitrary representation 
of the color group, and the transverse coordinates omitted, for simplicity), 



.11 V} 3 = g(T a V% , .11 Va = -g{VT% 3 , 



J a L 0(V, V ] ) = g{T a V r ) 



dO 



g{VT% 



dO 



(6.18) 



together with similar relations for the action of J R . These equations, as well as Eq. (6.10), 
can be rewritten as commutators between formal operators: 



J a L (x),Vl(y)] = g{T a V\x)) bc 5^\x 



y) 



Vl(y)] = g{V\x)T a ) b J 2 \x - y) , 
The left and right Lie derivatives span two independent SU(iV c ) algebras : 

J a L (x), J b L (y)} = -if abc J c L (x) 8W(x - y) , 
J a R (x), J b R (y)] = if abc J c R (x) 5^(x - y) , 



(6.19) 



J a L (x), J b R {y) 



0. 



(6.20) 



The first two equations above follow by enforcing the Jacobi identity in Eq. (6.19), whereas 
the last commutator, namely 



J a L (x), J b R (y) = J a L (x), V bc {y) J c L {y) + V bc {y) J a L (x), J c L (y) = 0, 



(6.21) 



vanishes because the infinitesimal gauge rotation of the Wilson line is compensated by 
the commutator of Jl with itself. 



We are now prepared to characterize the canonical Hamiltonian structure associated with 
the J1MWLK evolution: The Wilson line V bc (x) and the left Lie derivative J£(cc) play 
the roles of the canonical coordinate and its conjugate momentum, respectively, and the 

7 Note that the "left" and "right" nomenclature refers, by convention, to the action of the 
Lie derivatives on V* ; for the corresponding action on V, the "left" and "right" would be 
interchanged (see, e.g., Eq. (6.18) below). 
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Hamiltonian Hji MWLK [J l , V^] - as defined by Eqs. (6.16) and (6.17) — acts on the re- 
spective phase-space via the Poisson brackets defined by the 'left' part of the commutation 
relations 8 in Eqs. (6.19)-(6.20). The evolution of an arbitrary observable 0(V, V^) is then 
governed by the canonical equation of motion : 

(6.22) 

which after averaging with the target weight function (now expressed as a functional of 
V and V*) is equivalent to Eq. (3.8). 

In particular, when acting on dipole scattering operators so like Eq. (2.2), Eq. (6.22) 
produces the same evolution equations as originally obtained by Balitsky [10] from the 
analysis of projectile evolution in a strong background field. This is natural since the 
gluon mergings in the target wavefunction, cf. Fig. 2. a, can be reinterpreted as splittings 
in the projectile wavefunction followed by multiple scattering between the products of this 
splitting and the strong target field. In fact, the effective action (6.5) can be understood 
as describing one-step quantum evolution in the scattering between a high-density target 
(represented by the strong classical color field a) and a low-density projectile (the source 
of the weak field A~). With this interpretation, Eq. (6.5) is equivalent with the 'real' part 
of the effective action obtained by Balitsky (see Eq. (41) in the second paper of Ref. [9]). 



-00/0) 



H 



JIMWLK 



7 The general effective action 



Turning to the general case in which both background fields are strong (at least, for- 
mally), it turns out that the calculation is most conveniently performed by working in 
the background LC-gauge and using the expression (4.19) for the effective action. Indeed, 
the corresponding background field propagator is explicitly known (cf. Eq. (5.10)), and 
involves only one type of Wilson lines. Also, the vertex D~F +% = D~d + B l which appears 
in Eq. (4.19) is a total derivative (in a gauge-covariant sense) with respect to both x + and 
x~ , so the corresponding integrations can be easily done. In fact, the propagator (5.10) 
satisfies D~G l i (x,y) = and d+G i: >{x,y) = 0, so both integrations receive contributions 
from the boundary terms alone. 

After performing the integrations over x + and y + , one finds: 

x~,y-,x 

- F£(x~, x)W 00 ,_ 0C (x)F^ 00 (y-, x) - F^(x~, x)W^_ 00 {x)F+ i {y- , x)} . (7.1) 



8 Of course, the Poisson brackets of with itself, or with V, are trivial: [V^ b (x), V^ d (y)] = 0, 
etc. 
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By using F +l = d + B l and recalling that the field B l vanishes at x = — oo, one can also 
perform the integrations over x~ and y~, with the following result: 

X 

-- BUx)W^^x)B^(x) - B^(x)wL^(x)BUx)}. (7.2) 
where B l x+ (x) = B l (x + ,x~ = oo,x). 

As anticipated, the effective action (7.2) is completely determined by the asymptotic 
values of the fields at large x + and x~ , which in turn are the same as the respective 
fields along the edges of the interaction region in Fig. 3 (since the fields remain constant 
outside this region). This property is essential for the success of the renormalization group 
analysis, as it insures that the separation of scales is preserved by the quantum evolution: 
The effective dynamics at the soft scale 'sees' only a coarse-grained (in x + and x~) version 
of the dynamics at faster scales, but not also the short-range details of the latter. 

Our final step will be to express the action fully in terms of Wilson lines. To that aim, 
one should recall that the field B % is a pure gauge given by 

Bl + (x)= t -V x+ d*V^ (7.3) 

where = U^(x~ = oo) with the Wilson line given in Eq. (4.16). The longitudinal 
Wilson lines are explicitly known in terms of the field a in the Coulomb gauge, so it is 
convenient to fully rewrite Eq. (7.2) in this gauge. From Eq. (5.13), we deduce that 

VFoo,-oo(oo, x) = V^ (x)W O0 - 0O (x)V- O0 (x) 

= Pexp^igJ dz + 5A (z + , x~ = oo, cc)j , (7.4) 

whereas W 00) _ 00 (x) can be identified with the Coulomb gauge temporal line at x~~ = — oo: 
Woo-oofa) = W / oo,-oo( — oo, x). Let us introduce the simpler notations: 

W^oo,-oo(oo,a3) = Woo (a;), .^(-oo, x) = W^x). (7.5) 

By using Eqs. (7.3)-(7.5) and some simple manipulations, one can rewrite Eq. (7.2) in 
terms of the Wilson lines which delimitate the interaction region in Fig. 3 : 



AScfr = ^-JL^. Ini | Tr^X^Vl) + (d^WVl^) 

- (d'V^W^d'V^Wl - (d'V^W^d'vl^WlJ, (7.6) 



where we have also made use of the identity 



W ab = Tr[r a WT b W f ]. (7.7) 
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Note that the lower subscripts ±00 refer to x + in the case of V and V\ but to x" for 
W and (see also Fig. 3). Let us now discuss some general properties of the effective 
action (7.6), with emphasis on gauge symmetry and self-duality. 

The action (7.6) is invariant under ^-independent gauge transformations, which are the 
residual transformations permitted by the Coulomb gauge. Indeed, if h x +(x~) e SU(iV c ) is 
such a transformation, then W± oc — > h 00 (±oo)W± 00 h_ 00 (±oo) , — > h^— 00)140/1^(00) 
and — > /i 00 (oo)V ^/iJ )0 (— 00), so that the trace in (7.6) remains unchanged. 

From the discussion in Sect. 4 we expect the expression in Eq. (7.6) to be 'self-dual', 
i.e., invariant under the transformations (4.25). However, this symmetry is not manifest 
at the level of Eq. (7.6), where the temporal and longitudinal Wilson lines enter on a 
different footing. This asymmetry reflects our prior use of the background LC-gauge, in 
which the gluon propagator involves only the temporal Wilson line. However, after some 
integrations by parts, it is possible to rewrite Eq. (7.6) in such a way that self-duality 
becomes manifest. We shall write, e.g., 

Tr^yj^O = \j TrhO^KL - V^&V^)], 

X X 
X X 

and then repeatedly use Eqs. (7.4)-(7.5) to finally deduce two different but equivalent 9 
expressions for the effective action: 



AS, 



cff 



2ng 2 N c 



In 



Tr 



+ (VWtii&V&W-ooV-oo + Wiei&V&i&W-oJV-c 



(7.8) 



and respectively 



AS cS = 



27T 9 2 N C 



In 



Tr 



+ (d'W^WV^W^Vl^ + ^! oo (9V 0O )(9W 0O )y_ t J. (7.9) 



It is likely that this 'degeneracy' in the form of the effective action corresponds to the two 
original versions of it, Eqs. (4.19) and (4.20), respectively. In any case, it is now easy to 
check by inspection that Eqs. (7.8) and (7.9) get interchanged with each other under the 
following transformations 





v oo? 



H/_oo <— > VI 



t 



(7.10) 



9 The integrands in the two following expressions are Hermitian conjugate to each other and 
separately real, so they are indeed identical. 
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which represent a stronger version of the transformations (4.25) in which the ordering of 
the Wilson lines is also specified. Thus, the effective action is invariant under the duality 
transformation (or 'self-dual'), as anticipated. 

Eq. (7.8) has other interesting symmetries. For example, it is invariant under the following 
transformation 

W-oo-^-oo, V-o.-W*, W£-V£, Vl^W^, (7.11) 

which can be recognized as the clockwise 90° rotation of the diamond-shaped interaction 
area in Fig. 3. 

Note also that Eqs. (7.4)-(7.5) imply = V^W^V^, which can be rewritten as: 

VlW^V^Wl = 1 <=► J J- Y TrW <> = 1, (7.12) 

where W<> is the Wilson loop built around the interaction 'diamond' in Fig. 3. The fact 
that the overall Wilson loop is trivial is a consequence of the condition F + ~ = which 
expresses the separation of scales in the problem. (In fact, Eq. (7.12) is still another gauge- 
invariant expression of this separation of scales.) Eq. (7.12) can be also recognized as the 
generalization of Eq. (6.11) to the case where both background fields are strong. Indeed, 
as we shall shortly discuss, and W-oo play the same role in the general Hamiltonian as 
the functional Lie derivatives (6.9) in the JIMWLK Hamiltonian. Because of the condition 
F + ~ = 0, the (generalized) derivatives at the 'end points' x~ = oo and x~ = — oo are not 
independent of each other, but rather they are related by a gauge rotation, as manifest 
in either Eq. (6.11), or Eq. (7.12). 

By exploiting the self-duality property, one can deduce one more form for the effective 
action, which will be useful for comparison with the JIMWLK action in the dilute regime. 
Specifically, by applying the duality transformations (7.10) to Eq. (7.6), one finds 



x 



(9Wj o )(9W 0O ) + (9Wl 00 )(9W_ 00 ) 

- (9W 00 )Fi 00 (9Wl 00 )K - (dWl)Vi(dW^)vJ. (7.13) 



Any of the equivalent equations (7.6), (7.8), (7.9), and (7.13) defines a two-dimensional 
field theory in the space spanned by the asymptotic Wilson lines, with a local action. 
Because of the condition (7.12), only three among the four Wilson lines which enter any 
of these expressions should be treated as independent field variables; these can be chosen, 
e.g., as V£, V.oo, and W^. 

In order to complete the construction of the Hamiltonian theory, one still needs to specify 
the commutation relations among the independent Wilson lines. In principle, these could 
be deduced by using the expressions of the Wilson lines as path-ordered exponentials 
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of a a (x) and 5/5p a (x), respectively, together with the elementary commutation relations 
among these Lie-valued fields and derivatives; e.g., the Poisson equation (4.17) implies: 

= S ab S(x + - y + ) 5(x- - y-) A(x - y) , (7.14) 

with A(x — y) the two-dimensional Coulomb propagator introduced in Eq. (3.4). In fact, 
in the relevant Wilson lines, the field a a (x) appears only for asymptotic values of time, 
x + = ±oo; similarly, S/Sp a (x) enters only for x~ = ±oo (which truly means x~ and Xq in 
the notations of Eq. (6.9)). Accordingly, the commutation relations among and 
should receive non-trivial contributions only from their end points at x + = x~ = oo. 

But the construction of the general commutation relations is complicated by the fact that, 
e.g., the end-point derivatives S/Sa^ do not commute with each other, but rather obey the 
SU(iV c ) algebra in Eq. (6.20); via Eq. (4.17), this implies non-trivial commutators among 
the derivatives 5/5p£ Q . In the next section we shall discover that, in the dilute regime, 
the field variables should be treated as non-commuting variables (and similarly for 
a -oo) [38]. Hence, in the general case, the (end-point) gauge fields do not commute with 
themselves, neither do so the respective functional derivatives. If taken at face value, such 
non-commutativity properties at the level of the Lie-valued fields and derivatives would 
imply that a Wilson-line commutator like 

is not only non-vanishing, but also non-local in x~ ! This would spoil the two-dimensional 
character of the effective theory at high energy. We do not know yet what could be the 
solution to this difficulty in the general case (if any !), yet the fact that its two limiting 
cases — the high-density regime, cf. Sec. 6, and the low-density one, to be discussed 
in Sec. 8 — can be unambiguously formulated as two-dimensional, Hamiltonian, field 
theories strongly suggests that a similar construction should exist also in the general case. 
We postpone this issue for latter investigations. 

Let us finally check that, as anticipated at several places, the general effective action re- 
duces to the JIMWLK action in the high-density regime where the fluctuations become 
unimportant and the radiated field SA can be treated as weak. In this regime, the tem- 
poral Wilson lines and H 7 -^ can be expanded out in perturbation theory, whereas 
the time-dependence of the longitudinal Wilson lines can be ignored: Vl^ = V\ 

etc. To study this limit, it is convenient to start with the form (7.13) for the effective 
action, since there all the temporal Wilson lines appear under transverse gradients. Thus, 
to obtain the expansion of Eq. (7.13) to quadratic order in SA , it suffices to expand each 
temporal line to linear order: 

PW^ta -ig& A~ (oo,x), « ig dM~(-oo, x), (7.15) 

where A~(x~,x) = / dx + SA (x). Then, clearly, Eq. (7.13) reduces to Eq. (6.5), which is 
the effective action corresponding to the JIMWLK Hamiltonian. In the next section, we 



Sp a (x) 
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shall similarly consider the limit of the general action in the dilute regime, and construct 
the associated Hamiltonian theory. 



8 The low density regime: Bremsstrahlung Hamiltonian 



In the case of a dilute target, the color field a is weak and the longitudinal Wilson lines in 
the effective action can be expanded in perturbation theory. If one starts with Eq. (7.6), 
it is sufficient to consider this expansion to lowest order, e.g., 



t ao {x)^ig J dx-d*a a ±oc (x-,x)T a = ^^P±MT\ (8.1) 



z 



where we have also used the equation (4.17) relating the field a to the Coulomb-gauge 
color source p. In the last equality above, p(z) stands for the color charge density in- 
tegrated over the longitudinal coordinate x~ . Within the same approximations, one can 
neglect the difference between quantities in the Coulomb and the LC gauges, and also the 
£~ -dependence of the temporal Wilson lines (since the gauge rotations in equations like 
(7.4) become negligible for weak a); that is, VF_oo « = W. It is then straightforward 
to derive the following limiting form for the effective action [38, 39] : 



(2tt)* 



1*4 / ic xyz [pU^pM + P^Wp-ooiv) 



xyz 



pUxW^pl^-p^WW^pliy)], (8.2) 



with the transverse kernel JC xyz as defined in Eq. (6.7). The associated evolution Hamil- 
tonian -^brem is then obtained by replacing the color fields 5A~ in the temporal Wilson 
lines by functional derivatives with respect to the color source p a (x + ,x) = J dx~p a (x), 
cf. Eq. (3.14), which yields 

W(x) = Tex P (- 9 /d,+ ^A_), (8.3) 

The ensuing Hamiltonian has the structure anticipated in Eq. (3.15) and describes brems- 
strahlung in the BFKL evolution of a dilute target (cf. Fig. 2.b). As originally observed 
in Ref. [39], -£/brem is dual to the JIMWLK Hamiltonian (6.6) in the sense of the trans- 
formation (3.16). As explained in Sect. 4, this duality reflects the fact that the same 
diagram, namely Fig. 2.b, can describe bremsstrahlung with decreasing k + in the target 
wavefunction, or gluon merging with decreasing k~ in the wavefunction projectile. 

In fact, this duality refers to the complete Hamiltonian structure, including the Poisson 
brackets. Specifically, by using Eq. (8.3), one can deduce the commutation relations 
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p^x), W bc (y)} =g(T a W(x)) bc 5^(x - y) 
pl^x), W bc (y)] =g(W(x)T a ) bc 6M(x - y), 



(8.4) 



which are similar to Eq. (6.19) and show that the color charges act as infinitesimal 
gauge rotations of the Wilson line W at its end points (or, equivalently, as functional 
derivatives w.r.t. the field 8A~ at x + = ±00). This in turn implies that, within the 
Hamiltonian theory defined by I^brem, P^> an d P-00 m ust be treated as non-commuting 
variables, which provide two independent representations of the SU(iV c ) color algebra : 



PooW> Poo (J/J 



-if^YM^ix-y), 
0. 



(8.5) 



Like in the JIMWLK case, the commutators (8.5) are enforced by Eq. (8.4) together with 
the Jacobi identity and the following relation (cf. Eq. (4.8)) 



showing that the color charge densities at the 'end points' x + = 00 and x + 
not independent quantities, as they are related by a gauge rotation. 



(8.6) 



-00 are 



We recognize in the equations above the 'dual' version of the commutators (6.19)-(6.20) 
for the JIMWLK problem, with the precise duality transformation being now: 



pKx) , J a R (x) <— > pfL^ix) , V\x) — W(x) . 



(8.7) 



In particular, the relation (8.6) between and p°_ ^ is dual to the relation (6.17) between 
J£ and J|,, and they both emerge as particular limits of the more general relation (7.12). 
Still as in the JIMWLK problem, it becomes advantageous to use Eq. (8.6) to eliminate 
p a _ ^ from the problem and thus arrive at the following Hamiltonian 



if, 



BREM 



(27T) 



vr) 3 / 



1 + w x wl - w x wt - w z wi 



ab 



pU 



y), 



xyz 



which involves only two independent variables, p^ and W ab = W ba , and is manifestly dual 
to the JIMWLK Hamiltonian (6.16). The evolution equation for an arbitrary observable 
0(Poo, W) is then obtained as 



BREM i O 

which is the analog of Eq. (6.22) in the JIMWLK case. 



.9) 



By virtue of the duality property it is furthermore clear that the BREM Hamiltonian 
endowed with the above commutation relations acts on the temporal Wilson lines in the 
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same way as the JIMWLK Hamiltonian (6.16) acts on the longitudinal Wilson lines. 
Accordingly, the evolution equations generated by -£/brem for the operators built with W 
and (e.g., TrW (x)W jf (y)) are formally identical to the Balitsky equations [10] obeyed 
by the corresponding operators built with V and VL This should not come as a surprise: 
when focusing on the correlators of W, we are viewing the diagram in Fig. 2.b upside 
down, as the JIMWLK-like evolution of the projectile (our left mover). Then a 2-point 
function like TrW(x)W'(y) describes the scattering between a right-moving dipole and 
the left-moving system which is evolving. 

But, of course, our main interest when using Eq. (8.2) is not in describing the JIMWLK 
evolution of a dense projectile, but rather the fluctuations associated with bremsstrahlung 
in the evolution of a dilute target. These fluctuations can be probed via scattering with a 
projectile which is more complex than a simple dipole (since the projectile must measure a 
higher-point correlation function of the gluon fields in the target). Consider, for instance, a 
projectile made with two dipoles. According to Sect. 2, the amplitude for the simultaneous 
scattering of the two dipoles in the dilute regime is computed as: 

(T ( - 2 \x 1: y 1 -x 2: y 2 )) T ~ {(a a ( Xl ) - a a { Vl )f (a b (x 2 ) - a b (y 2 )) 2 ) r , (8.10) 

where a a (x) = / dx'a^x^, x) is the color field in the target at the time of scattering, and 
is related to p^x) via Eq. (3.4) : a a (x) = J d 2 y A(x — y) p^y). Note however a subtle 
point concerning the definition (8.10) of (T^) T : Since the fields a a (x) do not commute 
with themselves, the relative ordering of the individual dipole operators within Eq. (8.10) 
turns out to be important : (T(l)T(2)) r ^ (T(2)T(l)) r ! It remains as an interesting 
problem to clarify the physical meaning of this potential ambiguity, and construct the 
evolution equation for (T^) T according to Eq. (8.9). 

But such ordering ambiguities are expected to disappear, at least, in the large- N c limit, 
where the color dipole picture is valid [26-28] and the evolution equations for the dipolar 
scattering amplitudes have been recently constructed [18,33-35]. In this limit, and for 
sufficiently large energy, the evolution generated by .Hbrem should reduce to the equations 
in Refs. [33-35]. 

Furthermore, all such ambiguities trivially disappear in the BFKL approximation in which 
the fluctuations are completely neglected (meaning that the Hamiltonian is evaluated to 
second order in the functional derivatives). To deduce the BFKL Hamiltonian from the 
bremsstrahlung one, it is convenient to first rewrite the expression within the square 
brackets in Eq. (8.8) as (1 — W X W^)(1 — W Z W^). Then to the order of interest it is enough 
to expand the various Wilson lines to first order: 

/oo 
dx + 5A-(x + ,x). (8.11) 
-oo 

The BFKL Hamiltonian is finally obtained by using the identification (3.14) and neglecting 
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the time-dependence of p 



BFKL — 



J JC xyz f ace f bde p a (x) 



xyz 



Sp c (x) 5p c (z) 



5p d (z) 5p d (y) 



Ay)- 



[8.12) 



This is recognized as the expression of the BFKL Hamiltonian which emerges naturally in 
applications of the CGC formalism to the dipole picture [28, 33, 35]. It is furthermore dual 
to the corresponding expression obtained from the appropriate expansion of the JIMWLK 
Hamiltonian (6.16) [28,41]. 



9 Conclusions 



In this paper we have constructed the effective Hamiltonian describing the evolution 
of gluon correlations in a hadron wavefunction with increasing energy in the leading 
logarithmic approximation. In its essence, this Hamiltonian describes BFKL evolution in 
the presence of gluon recombination and bremsstrahlung, and thus of 'Pomeron loops'. It 
thus provides the appropriate framework to follow the evolution of a hadronic system all 
the way up in energy, from a dilute initial state up to a high density state characterized 
by gluon saturation and the formation of a color glass condensate. 

A priori, the effective Hamiltonian governs the evolution of a single hadron wavefunction 
(say, the target in a high-energy collision), but it can be also applied to scattering once a 
factorization scheme is known (as is the case, for instance, for a simple projectile built with 
dipoles). The Hamiltonian has the fundamental property of self-duality, which reflects 
boost invariance and suggests that the present formalism should naturally allow for a 
symmetric description of the evolution in the high-energy scattering. 

The Hamiltonian has been obtained within a renormalization group analysis which ex- 
ploits the separation of scales in rapidity and thus entails a coarse-graining in both the 
longitudinal and the temporal directions. The RG analysis is successful because the sep- 
aration of scales is preserved by the evolution: the interactions responsible for the evo- 
lution are localized within a space-time region which is small on the resolution scale of 
the quantum gluons which are integrated out in a single step of the evolution. Because of 
that, and of multiple scattering, the effective action describing one-step evolution involves 
the target fields only through Wilson lines. More precisely, it involves three independent 
group-valued, two-dimensional, field degrees of freedom, which are any three among the 
four Wilson lines forming a Wilson loop with unit trace around the interaction region. 
Two of these Wilson lines are path-ordered in the longitudinal direction and describe the 
recombination processes responsible for saturation and unitarity corrections. The other 
two are ordered in time and describe the bremsstrahlung processes responsible for gluon 
number fluctuations and the associated correlations. Both type of processes play an es- 
sential role in the evolution towards high gluon density with increasing energy. 
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The present construction encompasses and generalizes previous approaches in the litera- 
ture, like the effective action approaches in Refs. [6-9], the Balitsky-JIMWLK equations 
[10, 11, 14, 15], the recently derived evolution equations with Pomeron loops (as valid at 
large N c ) [18,33-36], and the effective Hamiltonian for the dilute regime by Kovner and 
Lublinsky [38,39]. In most cases, the correspondence with such previous results has been 
explicitly demonstrated in the appropriate limits. For instance, in the situation where one 
of the background fields is weak, our effective action reduces to that derived by Balitsky 
[9] for asymmetric scattering. It turn, this situation encompasses two different physi- 
cal regimes: the high-density regime (the radiated fields are weak) where the JIMWLK 
Hamiltonian [11, 15] has been recovered, and the dilute regime (the target fields are weak), 
where the BREM Hamiltonian [39] has been shown to emerge. 

But some other correspondences are still to be explored. Our general action, as given 
by any of the equations (7.6), (7.8), (7.9) or (7.13), does not look similar to the action 
obtained in Ref. [9] for the collision between two high-density hadrons. Also, the relation 
to the effective action for reggeized gluons by Lipatov and collaborators [7, 8] remains 
unclear to us. Finally, it would be interesting to extract the large- N c limit of our present 
results and thus verify whether this is equivalent to the Pomeron effective theory of Refs. 
[34,35], as a priori expected. 

Another important open problem refers to the understanding of the Hamiltonian structure 
of the effective theory in the general case. The Poisson brackets that would be most naively 
written down — by the straightforward generalization of the corresponding analysis in the 
two (JIMWLK and BREM) limiting cases — would spoil the two-dimensional nature of 
the effective action (in the sense of restoring the longitudinal and temporal coordinates), 
and thus are unacceptable. Other 'quantization' prescriptions, like the construction of a 
path integral in rapidity (so like in Refs. [17,42]), are currently under investigation. 
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